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SUMMARY 


The CMC fluid mechanics computer program system is being developed to 
transmit the theoretical evolution of finite element numerical solution 
methodology, applied to nonlinear field problems into a versatile computer 
code for comprehensive flowfield analysis. This report is Volume II of a 
three volume set and presents data deck procedures for the CMC three- 
dimensional Parabolic Navier-Stokes (PNS) algorithm. General data 
procedures are introduced, followed by detailed description of a juncture 
corner flow standard test case data deck. A complete listing of the data 
deck is given in Appendix A, followed by a detailed explanation of grid 
generation methodology in Appendix B. Subsequent appendices present 
descriptive tubulations of all commands and variables available to the 
user. These are in alphabetic order with cross-reference numbers which 
refer to storage addresses. Volume I of this report is referenced for 
details of the theoretical foundation, development of the finite element 
3DPNS algorithm, and discussion of results for the juncture corner test 
case. The CMC computer program structure and description are detailed in 
Volume III. 
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INTRODUCTION 


Input facilities for the CMC computer program are highly sophisticated 
and greatly simplify data deck preparation and modification. The program 
sequentially scans data deck card images and operates on command name data as 
encountered. Numerical data required for each command operation are input in 
free format directly following the command card. Command operations can cause 
vectors to be filled, initiate a series of solution operations or specify 
output formats and titles. Command card sequence is quite flexible and care 
has been taken to ensure that most operations which must be performed sequen- 
tially are specifiable under one command name. 

The CMC data deck is divided into seven sections for description. Exclu- 
sive of machine related job control statements, the seven sections consist of 
a FORTRAN MAIN program and accompanying subroutines (specific to external 
corner flows), namelist data, geometric description, output format specifica- 
tion, boundary and initial condition data and solution directing commands, see 
Figure 1. Each of these data and its subset is preceded by a command card 
image which directs a program activity, and upon completion returns control to 
the next command data card. The program operates in a dynamic storage mode 
and the function of MAIN (data deck section one) is to allocate sufficient 
storage for the IZ array which is internally sized as a function of the number 
of finite elements requested for a specific problem. The namelist section of 
the deck (section two) is used to specify scalar, integer and floating point 
data utilizing the FORTRAN namelist option. The data are read in namelists 
NAMEOQ and NAMEQ2, respectively, and stored in the arrays IARRAY and RARRAY. 

The geometric description section (data deck section three) contains data 
required to generate a finite element grid suitable for the solution. One- 
dimensional and rectangular two-dimensional discretizations are formed by 
specification of grid refinement along the coordinate axes. More complex 
boundary shapes (nonrectangular) necessitate data specification along all domai 
boundaries and certain interior subdomain boundaries. Generated data consist 
of grid point coordinates for each generated node and a node connection table 
which defines the finite element domains. 











CMC employs a highly adaptive output routine which allows for data 
specification of the scalar and array variables to be printed, scale factors to 
be applied to each variable and titling information to head each variable list. 
Each is specified under a command name in the output section (section four. 

Fig. 1) of the data deck. The program operates in nondimensional units and 
data specified scale factors are utilized to provide flexibility of print units. 

A reference length parameter (REFL in NAME02) is also available to scale input 
and output coordinate data to a problem reference length. In addition, command 
names are available for specification of problem identifying titles to be printed 
at various strategic locations i.e., at the beginning of each set of printed 
output. The fifth section of the data deck contains the required boundary and 
initial condition data which is specified at the solution nodes. The finite 
element method easily handles mixed boundary conditions, hence, both fixed and 
Neumann type are allowed. Parameter tables, which are a function of the marching 
direction (i.e., axial pressure gradient) are specified in section six. Linear 
interpolation of these tables determines the value of the parameter at each 
computational station. The intended solution algorithm is input in Section 
seven of the data deck. Variation in solution procedure is obtained through 
command card input of the various program functional modules. Nonlinear, 
repetitive solutions are readily handled by the Linkcall list which is interro- 
gated at each integration step. 

The sections which follow give examples of some of the computer program 
features which greatly simplify data specification and manipulation of computed 
results. These are followed by a detailed description of the data requirements 
peculiar to data deck major sections. The data deck description follows for the 
specific juncture flow case , including sample print and plot data. 


CMC DATA PROCEDURES 

Program flow is controlled by sequentially scanning the data deck in 80 
column increments. Subroutine BDINPT scans a card for a command name and 
initiates a program function for the given command. A list of the allowed 
acronyms is given in Appendix C. Eight columns are reserved for the command 
names beginning in card column 1. Beyond column 8, other data such as non- 
dimensional izing scalars or control parameters associated with a specific name 
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may be specified. If the function of a particular command name is to read and 
store a set of numerical or literal data, the cards containing the data must 
directly follow the command card. Control is returned to BDINPT by terminating 
a numerical data string with a "T" or a "blank card". Literal data are terminated 
with the word "DONE" beginning in card column 1. 

Most of the numerical and literal data, other than command data specified 
in the CMC input deck, may be input in free format. Data delimiters may be 
blanks or commas, thus allowing for esthetic and meaningful arrangement of 
numerical data. Exceptions are namelist data which utilizes the standard 
FORTRAN namelist option and certain special card types which combine literal 
and numerical data. Several features which greatly simplify sequential and 
repetitive data specification are available in free format. For example: 


Repetitive Numbers 

12. 5*7. 

T 


Fills Array 

12. 7. 7. 7. 7. 7. 



Repetitive Sequence 
(one per card only) 

2(5. 2. 4. 



Fills Array 

5. 2. 4. 5. 2. 4 



Repetitive Blanks (a) 

10. 12. 3*P 22. 

T 


Fills Array 

10. 12. A A A 22. 



Increment by a constant 

5*150 10 

T 


Fills Array 

10 60 110 160 210 


- 

Exponential Notation 

6. 10.E-2 14.E-4 

T 

• 

Fills Array 

6. .1 .0014 
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Solutions which proceed in marching integration fashion may be restarted 
at specific intervals by allocating file space and specifying the command 
SAVETAPE in the data deck prior to initiating integration. This command causes 
all pertinent arrays to be written on the specified file at each solution print 
interval. File space is minimized by specifying a small repeat number which 
rewinds the file following the specified number of writes. For example: 

SAVETAPE 93 T Save data in unit 9 

will cause the data arrays to be written on unit 9 at each print interval. 

On the fourth write, the file will be rewound before writing, thus obliterating 
the first three sets of data. A restart of the solution is accomplished by 
issuing command RESTART. In this case the second integer indicates the data 
file at which to start reading. For example: 

RESTART 91 T Restart at first set of data 

will cause the first set of array data on unit 9 to be read into the storage 
arrays. Caution must be observed not to specify a restart set of data which 
is beyond the last write indicated in the initial run. 

DATA DECK STRUCTURE 
Main Program ; (CMC Deck Section I) 

The function of MAIN (Fig. 2) is to dimensional ize and initialize the IZ, 
IARRAY and RARRAY arrays and call the control subprogram BDINPT . The IZ array 
resides in common block ARRAYS and contains all data array variables required 
for intersubroutine communication. The size of IZ must be overestimated 
initially since the arrays residing in it are dynamically set in the GO step 
as a function of specific problem size and solution type. On subsequent runs 
the size may be adjusted to fit the problem storage requirements. Actual 
storage required by the program is computed and stored (DECIMAL) at location 
100 in the IARRAY array. The IARRAY and RARRAY arrays reside in common block 
VARBLE and are utilized to conmunicate scalar integer and floating point data, 
respectively, between subroutines and I/O. All 5QC locations in each of 
these arrays are printed upon encountering command ICOND, thus permitting demand 
review of the many keys and parameters utilized throughout the program. 
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PROGRAM CORNER ( INPUT » OUTPUT * RSTRT » TAPE5= INPUT * 
* TAPE6=0UTPUT *TAPE3=RSTRT ) 

£ _ _ _ £ _ _ Jyj _ Q - - ^ 

COMMON / VARBLE / IARRAYC00500) » RARRAY(00500) 
EQUIVALENCE ( IARRAYC00092) * IZSIZE ) 

COMMON / ARRAYS / IZC 110000) 

COMMON LIST (200) 

NZSIZE = 110000 

CALL ERRSET < 207* 500* -1* 1* 0* 217 ) 

CALL ZEROTK 
100 CONTINUE 

CALL RESET < 00500* I ARRAY » 0 ) 

CALL RESET ( 00500 * RARRAY * 0*0 ) 

IZSIZE *> NZSIZE 

CALL RESET ( IZSIZE* IZ* 0 ) 

CALL RESET < 00200* LIST* 0 ) 

CALL BDINPT 
(30 TO 100 
END 


Figure 2. - CMC Main Program. 
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Namelist Control and Reference ( CMC. Deck Section II) 

All scalar data of consequence are stored in the arrays IARRAY and RARRAY 
(Appendix C). Many of these are FORTRAN namelist specifiable and are input 
under &NAME01 and &NAME02, respectively. Namelist read is initiated through 
command FENAME which initializes all storage locations to default values prior 
to reading namelist. The namelist tables also contain the name IARRAY and 
RARRAY, thus allowing data to be read into all 500 locations in each array 
regardless of namelist name specification. For example, NEQKNN = 5, can 
also be specified as IARRAY (58) = 5, under &NAME01I in Section II of the data 
deck. 

An alternative method of reading data into the scalar arrays is to use the 
IARRAY and RARRAY command names. These can be inserted anywhere in the data 
deck, external to namelist, and are useful for resetting print codes and dynam- 
ically calculating scalars to be used in subsequent steps. For example: 

RARRAY 23 1.5142 -3 T RESET TIME 

stores a new value of 1.5142 7 ALC into location RARRAY(23). The (-3) indicates 
division by ALC which is the variable stored in RARRAY(3). Likewise for integer 
data 

IARRAY 281 5048, 282 0, 283 71, 284 72 

stores 5048, 0, 71, and 72 in IARRAY locations 281 through 284, respectively. 

Note that these data cannot be continued on following cards. This does not restrict 
the method, however, since the next card could contain a similar command. 

Dynamic dimensioning is accomplished by filling the IZ array using the size 
specification N0DE(decimal ) as read under namelist NAME01. The command for this 
function is FEDIMN and it must be completed following NAMELIST read in Section 
II of the data deck. The dimensioning parameter NODE is discretization size 
oriented and is set slightly larger than the number of grid points in the 
solution. The starting addresses of the dimensioned variables in the IZ array 
are stored in the first NIZS (in NAME 01 ) locations in IZ and are printed in 
decimal form if KDUMP = 1 , in NAME01 . 
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Geometric Description (CMC Deck Section III) 

Finite element discretization of the flow domain is automatically generated 
from a coarse grid description. 

The parabolic-elliptic form of PNS requires a grid of space dimension one 
less than the solution geometry. The most complex grid description required, 
therefore, is two dimensional. A solution is obtained by "marching" the entire 
two-dimensional grid in a third space dimension with solution dependent step- 
size forming a partial adaptive grid. Boundary shape changes are controlled by 
geometric transformation parameters specified in data deck Section VI and are 
reflected in the solution through differential equation metrics. 

Simple Geometries 

A rectangular domain having grid lines lying parallel to the coordinate axes 
is the simplest geometry, and therefore, is simplest to input. Grid size vari- 
ation in each direction is accomplished using a geometric progression of finite 
element size in each direction. In the x 2 direction for example: 

m+1 . ? 

Y • + S Z P W (1) 

1+1 0 j=2 

where p is the specified geometric progression ratio and m is the number of 
finite elements to be generated in the x 2 direction and scaled by S. Since i 
in equation (1) is always positive, progression ratio greater than unity will 
form a grid which becomes coarser as x 2 increases and p equal to unity forms a 
uniformly spaced grid. The equation is applied piecewise over the elements 
allowing different progression ratios over each segment. For two-dimensional 
domains, equation (1) is applied simi lari ly in the second direction. 

The data description for rectangular discretizations proceeds by first 

subdividing the domain into smaller rectangles (subdomains) each of which 
may have a different progression ratio in each coordinate direction. The four 
subdomain case, illustrated in Figure 3a admits a fine discretization along 
the solid surface and in the tip region. Data for this case appear as follows 
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LINK2 


14 


VX1SCL 

-5.0, 

5 

0 . 

• 8, 

10 

1.0 

1.25 

VX2SCL 

-4.0, 

5 

0 . 

• 8, 

10 

1.0 

1.25 

NDECRD 

1 16, 

ELEM 

DONE 

1 

16, 

0 

T 




The first card is the discretization command card which effects a call to subr. 
DSCRTZ . The VX2SCL and VX1SCL commands read the subdomain data for the res- 
pective x 3 and x 2 directions of Figure 3(a). Each data string contains first 
coordinate, number of elements this subdomain, last coordinate this subdomain, 
and, progression ratio. Since the subdomain must be contiguous, the first 
coordinate of subsequent subdomains coincides with the last coordinate of the 
one preceding it and, therefore, is not repeated. The NDECRD command admits 
a choice of grid refinement to be used for solution, and generates the grid 
points. The data indicate first node in the x 2 direction (Fig. 3(a), last 
node in the x 2 direction, first node in the x 3 direction, and last node in the 
x 3 direction. In the example, 15 elements are described in each direction 
(VX1SCL) and (VX2SCL) and, hence, (1 16) will yield all elements described. 

If, however, only the upper half of the discretization were desired for a 
particular solution, the NDECRD data would appear as ( 6 16, 1, 16, 

0 T). The ELEM command causes the finite element node connection table to be 

formed. The number of elements and nodes generated (NELEM and NNODE) in IARRAY 
are determined and stored automatically in the discretization process. The 
number of generated nodes must be less than NODE specified in Namelist NAME01. 

Complex Geometries 

Secondary flow geometries having boundaries of general curvature and 
multiple discontinuities require a more sophisticated approach. The CMC 
multidimensional grid generator employs a local parametric coordinate trans- 
formation technique over independent quadrilateral and triangular segments of 
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the solution domain. A refined grid is generated over each local (single 
curved) subdomain and transformed into the global geometry. Figure 3(b) 
illustrates the method for a viscous corner flow problem where the potential 
core forms a curved boundary. Here four isoparametric quadrilateral subdomains 
having biquadratic variation are used to generate the two-dimensional triangu- 
lar mesh (see Appendix B). This method requires the specification of only 
coordinates in any order and a subdomain connection table. The connection 
table order is: corner grid points, side grid points counter-clockwise starting 

anywhere. A local coordinate convention is used to order the generated data. 

The origin is the first grid point in the connection table, the abscissa runs 
along side 1, between the first and second grid points in the table and the 
ordinate runs along side 4 between the first and fourth grid points in the table. 
The generated discretization is ordered in adjacent rows running between sides 
4 and 2 beginning at the local origin. For example, if subdomain 1 in Figure 
3(b) is numbered 

15 1 3 17 , 10 2 11 16 

and 5 elements per side were requested, the generated nodes would be numbered 
1 thru 6 along side one (15 1), 7 thru 12 along the adjacent row, etc. Note 
in Figure 3(b) that the number of generated elements along side (3b) must be 
the same for subdomains I and II in order to maintain nodal connectivity among 
the generated elements. 

The nonuniform grid exhibited in Figure 3(b) is obtainable utilizing two 
distinctly different approaches. The first method involves movement of the 
subdomain side grid points (Fig. 3(a)) toward the region of refinement. Since 
the number of generated nodes is equally distributed about the side node, a 
grid compression occurs for off-center side grid points resulting in a smoothly 
skewed grid over an entire subdomain. This method allows for great flexibility 
and control over spacial distributions of generated grid density with minimal 
input. The second method involves additional specification of geometric pro- 
gression ratios along subdomain sides and setting scalar variable JPR = 1 in 
the IARRAY array. This separates the grid refinement specifications from 
geometry input, thus providing a means for changing grid refinement without 
having to locate a new geometry point. It also overcomes, at little expense, 
the further drawbacks of method 1, involving nonuniqueness and boundary shape 
inaccuracies of highly skewed discretizations. 


11 


To use this option it is necessary to specify the number of subdomains 
(NSELEM) and grid points (NSNODE) in Namelist NAME01. Node must also be set 
in NAME01 larger than the number of nodes being generated. 

Report Format Output (CMC Deck Section IV) 

The output section of the CMC data deck is devoted to organizing the 
print formats. Several print command names are available, each having a 
different function. Table 1 lists the print commands and gives a brief des- 
cription of the function of each. Cards following the COMTITLE command are 
read and stored for print on command The command DESCRIPT 204 stores 

problem identifying titles for print at the beginning of each output header. 

The variable specified print formats in CMC allow the user to arbitrarily 
select and title the scalar and array variables to be printed. Scalars are 
normally grouped together in a header print which precedes print of the vari- 
able arrays. The titles are input under command DESCRIPT 332 and are formatted 
to fit 5 titles per card with 16 characters maximum per title. The scalar 

variables associated with each title are input under command IONUMB. These 

data consist of integers indicating the locations in RARRAY to be printed. 
Multipliers to be applied to each printed number are likewise input as integer 
RARRAY locations following the command MPARA. The example below illustrates 
the input which prints velocity in the English, MKS, and CGS systems of units: 

DESCRIPT 332 

VELOCITY FT/SEC METERS/SEC CENTIMETERS/SEC 

DONE 

IONUM -1 

200, 3*27 T 

MPARA -1 

2*2, 164, 163 T 

The DESCRIPT 332 command stores the four titles in the order read. The vari- 
able to be printed under each title, in this case RARRAY (200) = 0.0 and RARRAY 
(27) = , are listed under the IONUM command and zeros are printed as blanks. 
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Table 1 Print Commands 


Command 

Function 

COMTITLE 

Reads title which is printed below the CMC ! 
symbol 

DONE 

Terminates literal data 

DESCRIPT 204 

Solution print heading 

DONE 

Same as above 

DESCRIPT 332 

Solution print parameter titles 

MPARA -1 

RARRAY locations of multiplier 
TO BE APPLIED TO PARAMETER PRINT 

IONUMB -1 

RARRAY locations of parameters 
to be printed 

DESCRIPT 203 

Titles to head dependent variable 
print 

IOSAVE -1 

List of dependent variables to be 
printed. 

IOMULT -1 

List of locations in RARRAY of multipliers 
to be applied to each dependent variable 
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The printed scalars are scaled by data in RARRAY locations specified under the 
MPARA command. In the example, therefore, U m is printed three times, scaled 
by RARRAY(2) for print under heading FT/SEC, scaled by RARRAY(1|64) for print 
under heading METERS/SEC, etc., where the multipliers are 1.0, feet to meters 
and feet to centimeters conversion factors. A typical header print obtained 
upon command LINK2 6 is illustrated in Figure 4, together with data descrip- 
tion which created it. Note that the (x-x-) in Figure 4 print is specified as 
999 under command IONUMB. 

Titles heading the print of array variables are input under command DESCRIPT 
203. Input format is 16 columns per title, 5 titles per card; 1 title is 
input for each variable to be printed. The arrays to be printed under each 
heading are assumed NNODE long and follow command IOSAVE. The program array 
variable/names are stored in the first NIZS locations of the IZ array and are 
described in the Section V data deck description following. Multipliers used 
to dimensional ize the print are specified for each variable following command 
IOMULT. Integer specification of RARRAY locations is used as described above 
for scalaHdata print. 

Print format control is obtained for array printing through namelist data 
specification. The maximum number print digits are specified in IARRAY(20) = 
NPRNT in NAME 01 . Output field size is controlled by specifying the required 
digit accuracy in IARRAY(22). Output for a typical case having simple geometry 
is illustrated in Figure 5. The format assumes leading decimal point mantissa 
(not printed) and the exponent is printed once for each variable directly 
following the array name. 

The print form reflects problem geometry and approximate coordinates are 
tabulated on the abscissa and ordinate of a Cartesian rectangular system for 
each variable printed. Coordinate values are in reference units and are real 
numbers with assumed leading decimal as described above. 
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DESCRIPT 332 T IOPAR PARAMETER TITLES FOR OUTPUT. 

REFERENCE ENGLISH-FT ENGLISH-IN M-K-S C-G-S 


LENGTH 

LJFI nr TTY 

.FT 

.FT/S. ...... 

• IN * 

» M 

A M/S 

.CM. . . 
.CM/S. 


DENSITY 

•LBM/FT3* ... 

*N» A. ....... 

KG/M3 

G/CC.. 

• ♦ ♦ ♦ 

TEMPERATURE . . . ♦ 

.RANKINE.... 

• N* A. ....... 

.KELVIN 

•N.A. . 


ENTHALPY 

♦ BTU/LBM. ... 

»N* A. ....... 

• KJ/KG 

. N . A . . 

...... 

FROZ. SPEC. HEAT. 

• BTU/LBM-R • • 


.KJ/KG-K.... 

.N.A. . 


VISCOSITY 

• LBM/FT-S. • . 

. N . A . ....... 

.NT-S/M2. ... 

♦ POISE 

LOCAL PRESSURE. 

. PSF 

. PSI ........ 

'KNT/M2 

♦TORR. 


LOCAL SOLUTION. 

TODCT Y . * 

.MACH NO.... 

. DPDX1 

. G22 ....... 

. G32 ........ 

.EPSILON.... 

.ENERGY 

.G23. ....... 

.INT. 
.FI * • • • 

VAR. .. 

kJi.mrnM u'Q- - . . * 

. H3 1 

. G33 



Xl/LREF 

• DX1/LREF. . . 

. DX1M/LREF. . 

• REFL 

RE NO. 


DONE 

MPARA -1 

5*2 , 2 2 162 164 163, 3*2 164 163, 3*2 170 174, 

3*2 165 2, 2 -175 3*2, 3*2 176 2, 3*2 177 178, 

2 2 169 160 167, 3*2 108 2, 5*2, 5*2, 5*2 T 

IONUMB -1 

999, 5*200, 999, 200 4*43, 200 27 200 27 27, 

200 10 200 10 10, 200 58 200 58 200, 

200 97 200 97 200, 200 30 200 30 200, 

200 38 200 38 38, 999, 39 4*36, 300 154 100 135 122, 
398 4*11 186, 200 4*11 139, 11 12 14 R5 47 T 


4(a). - Header Print Format Data, 


REFERENCE ENGLISH-FT ENGLISH-IN M-K-S C-G-S 


-x-x-x-x 

X 

i 

X 

1 

X 

1 

X 

1 

-x-x-x-x- 

x-x-x-x 

1 

X 

l 

X 

1 

X 

l 

i rur.TU , . . 

t FT r 

. IN 

, M 

.CM 


1 OOOOOOE+Ot 

1200000E+02 

3048000E+00 

3048000E+02 

uri nr tty . . . _ . 

* FT /S 

N ft | , x l J; 

.M/S 

.CM/S 


1000010E+03 


3048000E+02 

3048000E+04 

DENSITY 

.LBM/FT3. . . . 

.N.A. ....... 

KG/M3 

G/CC 


2330814E-02 


3733965E-01 

3733965E-04 

TEMPERATURE. . , . 

♦RANKINE.... 

. N.A. »*»«... 

.KELVIN 

♦N.A. ....... 


5300000E+03 


2944444E+03 


ENTHALPY 

.BTU/LBM. ... 

•N.A........ 

.KJ/KG...... 

•N.A*. ...... 


1760684E+04 


4092533E+04 


FROZ. SPEC. HE AT. 

.BTU/LBM-R. . 

. N.A. ....... 

.KJ/KG-K, . . . 

.N.A. ....... 


7721760E+01 


3230784E+02 


UTQrntJTTY * * 

* i rm/ft-s * * * 

f N t A , 

♦NT-S/M2. . . ♦ 

.POISE 


4074412E-06 


6062725E-06 

6062725E-05 

[-X-X-X-X 

-X-X-X-X 

-x-x-x-x- 

X-X-X-X 

- X - X - X - 

LOCAL PRESSURE. 

.PSF 

.PSI. 

KNT/H2 

.TORR 

9964654E+00 

21 16800E+04 

1 469906E+02 

1013524E+06 

7601429E+03 

LOCAL SOLUTION. 

.MACH NO.... 

.DPDX1 

.ENERGY 

.INT. VAR... 

0 E+00 

8861539E-01 

0 E+00 

5847034E+00 

0 E+00 

T2PFIY . .. , , , 

. H21 ........ 

.fi?? 

, G23 

.FI .......... 

0 E+00 

1000000E+01 

0 E+00 

0 E+00 

0 E+00 

tjunrnM u'Q 

. H31 ........ 

t n32 

,G33 

,G1 


1000000E+01 

0 E+00 

0 E+00 

2129497E-1 4 

Xl/LREF 

.DX1/LREF . • . 

.EPSILON, . . . 

♦DX1M/LREF • • 

.REFL RE NO. 

1100000E-01 

1000000E-02 

1000000E-01 

0 E+00 

5720616E+06 

-x-x-x-x 

-X-X-X-X 

-X-X-X-X- 

X 

i 

X 

1 

X 

1 

X 

-X-X-X-X 


x-x-x-x-x-x-x-x 

N- N+ PASSES PRINT 

0 0 1 0 


4(b). - Header Print. 


Figure 4. - Typical Solution Header Data and Print. 
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In addition to the flexible format print described above, other less 
flexible print may be obtained which provides input validation and verifica- 
tion of problem definition. Setting KDUMP = 1 in NAME01 will effect a line by 
line data deck reflection, including each command card and its associated data 
as read. Figure 6 illustrates typical data verification print. Common data 
errors such as not specifying "T" or blank card delimiters are easily detected 
since the reader will continue to read cards, including command names, until 
it encounters one of these. The print subsequently appears as continuous data 
reflection of all the cards under the original command name. 

Scalar data stored in the arrays RARRAY and I ARRAY are; printed by inserting 
the command ICOND. Figure 7 illustrates the print format. This print provides 
verification of the namelist data specification, and program computed scalars 
which are stored in these arrays. 

Upon call to QKINT, a series of prints is initiated. The first, illustrated 
in Figure 8, prints the variable numbers and types in the solution. If a re- 
start unit is requested, the unit number is printed. The order of calls for 
solution process is listed and identified, followed by the variables and mul- 
tipliers to appear following the header page at each DELP solution print inter- 
val. Following this print, a call to ICOND is initiated to print the RARRAY 
and IARRAY scalar lists. This is followed by a node map print which matches 
node numbers with their coordinates. A subsequent standard solution print, as 
described above, is output which lists the specified and default initial 
conditions. 

Description of various debug print which is useful for verifying algorithm 
correctness and debugging the program is given in Volume III, the Programmer's 
guide. 

Boundary and Initial Condition Specification (CMC Deck Section V) 

CMC admits the simultaneous solution of multiple differential equations. 
Dependent and parameter array variables are stored sequentially, beginning at 
the address stored in IZ ( 48) and are input in arbitrary order following command 
IP I NT. The number of differential equations, hence, number of dependent vari- 
ables in a solution is specified by NEQKNN in NAME01, and the total number of 
variables specified in IPINT is input in NEQ. Each specified variable is 
actually a composite number which is decoded by the program for identification. 
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Figure 6. - Sample Data Reflection Print (KDUMP = 1 in NAME01) 


STATUS OF X ARRAY VALUES. 


1 

NO 

S 

1 

2 

MLTDIF 

= 

0 

3 

I ROW 

a 

1 

4 

NEYHTD 

a 

1 

5 

NOMEGD 

a 

0 

6 

KODCi 

a 

0 

7 

K0D5 

= 

2 

8 

KPRINT 

a 

0 

9 

NRSTRT 

a 

0 

10 

NO 

a 

0 

11 

NPTUOF 

B 

1 

12 

NCC 

= 

0 

13 

NCOORD 

a 

0 

14 

NFLEH 

B 

325 

15 

NN 

a 

0 

16 

NNDDE 

= 

190 

17 

NNS 

a 

0 

IB 

NPART 

a 

2 

19 

NGE 

a 

171 

20 

NPRNT 

a 

132 

21 

NROU 

= 

2 

22 

NC 

a 

8 

23 

ND 

a 

4 

24 

NPTEL 

a 

3 

25 

KIND 

a 

4 

26 

KOUNT 

a 

1 

27 

NSK1F 

a 

193 

28 

1PASS 

a 

0 

29 

1RUN 

a 

0 

30 

NF‘ 

= 

0 

31 

NEQ 

a 

8 

32 

IWt.SEP 

- 

19 

33 

NBUG 

a 

0 

34 

LPRINT 

* 

100 

35 

npsicc 

a 

0 

36 

NPRSCC 

E 

0 

37 

NEWPRT 

a 

0 

38 

NOUTVC 

a 

10 

39 

NlJRVDX 

a 

0 

40 

NOUTS 

a 

1 

*1 

N2UAKE 

a 

0 

42 

INPUT 

a 

5 

43 

NEQADD 

a 


44 

N ST ORE 

= 

0 

45 

NEXP 

= 

19 

46 

NF 

3 

4 

47 

LG 

a 

19 

48 

NDOF 

a 

0 

49 

NOEL 

Ti 

3 

50 

LCQL 

a 

30 

51 

NDBL 

= 

0 

52 

KROU 

a 

46 

53 

NH2 

a 

0 

54 

NHMALF 

a 

0 

55 

NODE 

a 

193 

56 

NDOUTZ 

3 

0 

57 

IASVEC 

a 

0 

58 

NEQKNN 

a 

5 

59 

NCPTAB 

a 

1 

60 

NMBOUT 

a 

13 

61 

KDUNP 

S 

0 

62 

NTITL 

a 

3 

63 

NIND 

a 

1047 

64 

NSh 

a 

10 

65 

NJ 

a 

0 

66 

NRHOSS 

3 

0 

67 

NS 

a 

26 

68 

NI 

* 

4633 

69 

NBSET 

a 

0 

70 

NDTSET 

a 

0 

71 

ISTATN 

3 

0 

72 

IVALUE 

a 

0 

73 

ISTAD 

a 

0 

74 

IRRAY 

a 

0 

75 


a 

0 

76 


3 

0 

77 


= 

0 

78 


a 

0 

79 

NQQ 

a 

9 

80 

NPRDCT 

*L 

0 

81 

IVUVAL 

= 

0 

82 

IVSTA 

a 

0 

83 

IFSLT 

= 

0 

84 

NVARD 

a 

172 

85 

NVAR1 

~ 

171 

86 

KENT 

K 

1 

07 

NRTAPE 

a 

0 

88 

NLINE 

* 

60 

89 

NEND 

a 

349 

90 

NY Y 

a 

4 

91 

NZZ 

3 

4 

92 

I2SIZE 

a 

*0000 

93 

NTK 

a 

2 

94 

NI2S 

a 

250 

95 

NVELTB 

a 

0 

96 

I BASE 

a 

200 

97 

ITKE 

a 

1 

98 

IHAT 

a 

1 

99 

IBL 

a 

1 

100 

IREND 

* 

66216 

101 

NPRESS 

a 

0 

102 

KOUT 

a 

0 

103 

MTM 

a 

0 

104 

NFL IP 

a 

1 

105 

IPTSF'L 

aa 

0 

106 

KK 

a 

19 

107 

NE1E2 

a 

0 

108 

NCNADH 

a 

8 

109 

NTAPER 

a 

0 

no 

NPRESH 

a 

0 

111 

I MAX 

* 

9 

112 

NWALLS 

a 

4 

113 

KPLVAR 

a 

0 

114 

HSSD 

a 

0 

115 

IFORCE 

a 

6904? 

116 

IG1NIT 

a 

0 

117 

IBUG1 

a 

0 

118 

IBUG2 

a 

0 

119 

NGDUG 

a 

0 

120 

IFR 

a 

0 

121 

NSPEC 

3 

4 

122 

I WRIT 

a 

0 

123 

IGAS 

a 

0 

124 

NDERIV 

a 

2 

125 

NCALLS 

a 

2 

126 

I SPEED 

a 

0 

127 

IDIFRT 

a 

-2 

128 

ICNTND 

a 

190 

129 

IMATAB 

a 

2 

130 

NGETH 

e 

1 

131 

NBQRD 

a 

52 

132 

IPWRIT 

a 

-12 

133 

I TOP 

e 

4 

134 

NPSIfeD 

a 

2 

135 

NPSBDL 

a 

0 

136 

KSAV 

s 

0 

137 

LPSIA1 

a 

0 

138 

LPSIAM 

a 

0 

139 

INGAS 

a 

0 

140 

NTCNTS 

a 

0 

141 

KFSL 

a 

0 

142 

NOUTPR 

a 

68 

143 

I CALLS 

B 

0 

144 

ITRAN 

a 

0 

145 

IMIN 

a 

0 

146 

NSFOBE 

a 

1 

147 

NPAST 

a 

0 

148 

NONC 

a 

0 

149 

IT 

a 

0 

150 

NODNO 

a 

0 

151 

KUFLXS 

a 

0 

152 

INITKE 

a 

0 

153 

NPUNCH 

» 

0 

154 

NSDFCF 

a 

0 

155 

NS2DFC 

B 

1 

156 

ITWOCL 

a 

0 

157 

NPASS 

a 

0 

158 

INITCN 

B 

0 

159 

IL 

a 

1 

160 

IH 

a 

1 

161 

NPVSX 

a 

2660 

162 

IEF’SET 

a 

0 

163 

NPSIST 

a 

0 

164 

NPSIND 

a 

0 

165 


a 

0 

166 

N3DPNS 

s. 

0 

167 

KNTF’AS 

a 

200 

168 

NDPRES 

B 

0 

169 

KCDC 

a 

0 

170 

NDC 

a 

0 

171 

NBCND 

a 

194 

172 

LOUD 

a 

1 

173 

NCOhOC 

= 

0 

174 

NCOMTD 

a 

3 

175 

IFSL 

a 

0 

176 

NTKS 

= 

10 

177 

NU2P0S 

a 

3 

178 

NU3P0S 

a 

3 

179 

LOC 

a 

1 

180 

NSTRT 

s 

11 

1Q1 

NVELP 

3 

0 

182 

LOGS 

a 

0 

183 

IRSLHS 

a 

0 

184 

NGTPRS 

a 

11 

185 

IP 

a 

0 

186 

I T DA 

a 

0 

187 

ITDB 

a 

0 

168 

NBUF 

a 

0 

189 

NSTIJ 

a 

8 

190 

NHGUT 

a 

2 

191 

NM 

3 

3 

192 

N2M 

a 

6 

193 

NM2 

a 

9 

194 

NDP 

= 

10 

195 

KFXBND 

a 

0 

196 

ITWALL 

a 

0 

197 

NTPRNT 

as 

2 

198 

NPGRDT 

a 

4 

199 

NPGRDV 

a 

4 

200 

INTPSI 

B 

0 

201 

NSTRFP 

a 

1 

202 

NSTRF 

a 

14 

203 

NTCRDM 

a 

1 

204 

INTOHG 

:= 

0 

205 

ITIhER 

a 

0 

206 

NhDL 

a 

8 

207 

IBOT 

a 

3 

208 

NVRHS 

a 

8 

209 

NBIPRS 

a 

0 

210 

NCOhPG 

a 

35 

211 

IBIDRV 

a 

0 

212 

HLTRHS 

a 

1 

213 

IPHIPR 

a 

0 

214 

LhLT 

a 

0 

215 

1SUPRS 

a 

0 

216 

ITPSCL 

a 

0 

217 

NOPSCL 

a 

0 

218 

NCLPSI 

a 

7 

219 

KR 

a 

0 

220 

NBCNDT 

a 

0 

221 

NCYL 

= 

0 

222 


a 

0 

223 

MLTLHS 

a 

3 

224 

ILHS 

a 

0 

225 

NORM 

a 

7 

226 


a 

0 

227 


a 

0 

228 


a 

0 

229 


a 

0 

230 


B 

0 

231 

IADSET 

s 

0 

232 

NPHOD 

a 

200 

233 

NODES 

a 

170 

234 

KNTARF 

a 

0 

235 

KNTSHP 

a 

0 

236 

N5AF 

a 

0 

237 

NARFPR 

a 

0 

238 

NSHAPF 

a 

0 

239 

NARF 

a 

0 

240 

NGRDSH 

a 

0 

241 

NWAKF 

a 

0 

242 

INAF 

a 

0 

243 

NODFT 

K 

0 

244 

NOD TOT 

a 

0 

245 

NAF 

a 

0 

246 

NLC 

“ 

0 

247 

NEC 

= 

1 

24B 

KFIX 

a 

0 

249 

JSHP 

a 

0 

250 

IPHI 

« 

0 

251 

NSNOOE 

a 

19 

252 

NSELEM 

a 

2 

253 

JCOORD 

a 

0 

254 

IPIN 

a 

0 

255 

IDIAGI 

a 

0 

256 

I PC 

a 

0 

257 

NFLUX 

a 

0 

258 

NFX 

= 

0 

259 

IS IDE 

= 

4 

260 

NVAR 

a 

3 

261 

IBC 

a 

0 

262 

NDCNT 

a 

0 

2*3 

LOCK 

= 

190 

264 

L.NEL 

a 

190 

265 

KSIDE 

a 

32 

266 

JSIOE 

= 

0 

267 

KTE 

a 

0 

268 

NDIUHF 

a 

0 

269 

NPOT 

a 

t 

270 

ITERPF 

a 

0 

271 

IPRINT 

a 

0 

272 

KSKIP 

= 

0 

273 

NSTAG 

- 

0 

274 

NU 

a 

0 

275 

NL 

a 

0 

276 

NUS 

a 

0 

277 

NPTS 

a 

101 

273 

NNPT 

-• 

0 

279 

NDIM 

a 

200 

280 

N5EL 

a 

2 

281 

IUYY 

3 

0 

282 

NV 

a 

0 

283 

N01 

a 

0 

284 

N02 

a 

0 

285 

NACA4D 

a 

0 

286 

NEQAV2 

=* 

1 

287 

NEGA03 

a 

1 

288 

NSORC 

- 

0 

2B9 

IPLOT 

a 

0 

290 

JQAD 

a 

0 

291 

JPR 

= 

1 

292 

JCYL 

a 

0 

293 

IAXSYN 

a 

0 

294 

IRAT 

a 

0 

295 

IPLOTV 

a 

0 

296 

IPCOEF 

a 

0 

297 

IPPSET 

a 

0 

298 

NPSETS 

a 

0 

299 

IPCFIT 

a 

0 

300 


a 

0 

301 

NIHPLT 

a 

1 

302 

NHCNTR 

a 

0 

303 

NITER 

a 

0 

304 

NELPAS 

a 

0 

305 

NBAND 

a 

23 

306 

NhBJAC 

= 

4370 

307 

NCONV 

a 

0 

308 

N3MR0U 

a 

386 

309 

N4NRGU 

a 

579 

310 

NFHROW 

a 

386 

311 

JORDER 

a 

1 

312 

IUONLY 

a 

4 

313 

NR 

E 

1 

314 

KQDE 

a 

0 

315 

NPT 

a 

0 

316 

NIT 

a 

1 

317 

NDBGPT 

a 

3 

318 

KKSKIP 

a 

0 

319 

K DP ASS 

a 

0 

320 

ICHI 

a 

0 

321 

IPRNT 

a 

0 

322 

ISTART 

a 

0 

323 

NLAST 

a 

0 

324 

NCNTIT 

a 

0 

325 

NOUEDG 

a 

0 

326 


a 

0 

327 

IDOXST 

a 

2 

328 

NVISC 

a 

136 

329 

NRJACB 

B 

1 

330 

LMDJAC 

a 

0 

331 

NCONSV 

a 

0 

332 

JACSAV 

a 

0 

333 

INOCRN 

a 

1 

334 

1TRMAX 

a 

0 

335 

MAXPAS 

a 

0 

336 

NPDDUG 

a 

40 

337 

NTRLST 

a 

0 

338 


a 

0 

339 


a 

0 

340 


a 

9 

341 


3 

0 

342 


a 

0 

343 

N0BQ23 

a 

0 

344 

IBL AS 

a 

0 

345 

IU1C0N 

a 

0 

346 

NSPRED 

a 

0 

347 


a 

1 

348 


a 

1 

349 


a 

0 

350 


a 

0 


Figure 7. - Scalar Arrays Print 
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3 VARIABLES BEING INTEGRATED * 

1 5 A 2 3 

8 VARIABLES IN SOLUTION. 

1562 3 78? 

ORDER OF CALLS AT END OF QKNUIN 

LINKLM 4> - CONTES 

PRINTOUT VARIABLES 

1248 2240 3248 1247 3248 6240 7248 3271 4271 

6271 7271 8271 

PRINTOUT VARIABLE MULT IPL IERS . 

2 2 2 21 2 2 2 -2 2 

-2 2 -2 

PRINTOUT VARIABLE FACTORS. <N> 


N 

.EQ. 

U 

STRAIGHT PRINT. 



N 

• EQ. 

21 

MULT. LAST VEC. 

BY THIS ONE 

- STORE IN LAST 

N 

.EQ. 

3 > 

ADD TO LAST VEC « 

- STORE IN 

LAST VFC. LOC. 

N 

.GE. 

4 i 

RAISE ENTRIES TO 

< N-2 ) POWER « 

N 

.LT. 

0> 

TAKE NTH ROOT OF 

ENTRIES IF 

.GT. 0.0 


1 


1 l 1 

1 

1 1 1 


1 


1 




5271 


2 


Figure 8. - Solution Sequence and Variable Definition Print 


The composite consists of two parts: 

0 000 

The first digit string contains the variable number and is specified 
only if the array appears as a variable under command I PINT . 

The second digit string contains the indirect address of the variable. 

The indirect address is the position in the address array which points 
to the variable of concern plus, the integer constant IZBASE stored in 
the IARRAY array. 

For example, specification of variable 1248 would cause integration of depen- 
dent variable No. 1 beginning at the address stored in I Z( 200 + 48) where 
IZBASE is 200. Names and brief descriptions of the program variables are 
given on page 32 and in Appendix E. 

Boundary conditions are separated into two types for clarity of description. 
The simplest to apply is fixed values of each dependent variable. The nodes 
where the variable is to be fixed are specified following the command KBN0 I, 
where I is input beyond column 8 and is the number of the dependent variable 
to which the boundary condition will apply. The boundary nodes may be speci- 
fied as literal data, therefore, if actual node numbers are to be input they 
must follow the literal data command ADD. For example, if nodes 1-5 are to 
be fixed for dependent variable 5 the input would appear as: 

KB NO 5 

ADD 

1, 2, 3, 4, 5 T 

DONE 

Note that a DONE command must be used to terminate literal data. For rectangu- 
lar discretizations, the input is simplified by requiring only one name to 
specify an entire side. The names BOTTOM, RIGHT, TOP, and LEFT identify the 
side to be fixed. Interpretation of the names is associated with their appear- 
ance for a geometry in the first quadrant of a 2D rectangular Cartesian coor- 
dinate system. Input for the previous example would appear as follows if nodes 
1 thru 5 represented the entire discretization along the abscissa. 
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KBNO 


BOTTOM 

DONE 


Solution domain gradient boundary conditions are specified at grid points 
as a a and a 3 coefficients as they appear in Volume I (Eq. 35). Data are 
specified under the KBNO command as described above with special format to 
allow for specification of the coefficients. On the KBNO card, a 1 follow- 
ing the dependent variable number indicates that gradient boundary conditions 
are being applied. The literal node specification cards, such as BOTTOM 
contains the specified ai and a 3 coefficients. For example; 

KBNO 1 1 

BOTTOM 0 1 5, -.1 2, 0. 2 

DONE 

applies an a^ coefficient of -.1 times RARRAY(.2) and an a 3 of 0. times 
RARRAY(.2) to the nodes along the abscissa of a rectangular domain beginning 
at the left boundary (0 displacement) and spanning the first five nodes. For 
the case of viscous flow boundaries near a wall, an internally generated shear 
stress is applied by specifying a 3 = 1.234. By specifying a = 4.321 a second 
internally generated gradient boundary condition is applied which essentially 
satisfies the global continuity equation. Thus, for a two-dimensional flow 
problem, a du 2 boundary condition is applied by simply letting a 3 = dui 
3T W 

Primary Flow Parameter Tables (.COMOC Deck Section VI) 

The Parabolic Navier-Stokes Equation system is marched in one direction 
and hence the grid in the marching direction is determined by the solution 
step size. 
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Since the step size is bariable and completely unpredictable, a linear inter- 
polation scheme is used to evaluate parameters which are a function of the 
marching direction. As an example, external flows over curved boundaries 
generate an axial pressure gradient in the flow direction. The pressure varia- 
tion may be solved for using an inviscid analysis and applied to a boundary 
layer solution through tabular input. Since linear interpolation is utilized, 
the table must contain sufficient data to effectively describe the pressure 
curve. The abscissa locations are specified in free format under command VX3ST 
and pressure levels are specified under command VPVSX. 

0. .05 .1 .15 .3 .5 T 

VPVSX 

1. -.6 -.5 -.3 -.1 -.05 T 

The above data complete a table which is interrogated for pressure level upon 
command LINK1 4. Note that NPVSX must be set to at least the table size in 
Namelist NAME 01 . (Vol . I, P. 16). 

Piecewise linear tabulated coordinate data describing variable geometry in 
the flow direction (Vol. I, P. 16) is input under commands VU3P0S (xj direction 
coordinates) and VU3VAL (F 21 (x x ) and F 22 (xj ) coordinates). The tabular 

functions are input in the order F 22 , F 21 as illustrated below. 

VU3P0S -1 

1.0 1.05 1.1 1.2 2.0 T 

VU3VAL -1 



Table storage arrays for variable geometry coordinates are dimensioned by 
setting NU2P0S and NU3P0S in namelist NAME01 equal to the number of table 
values input (i.e., for the above case =5). Note that variable geometry 
requires a nonzero ALC which may either be specified in namelist NAME02 or 
set internally by calling subroutine DIMEN (LINK3 4) prior to specifying 
variable geometry data. 
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Table look-up for the above is strictly interpolator and, therefore, 
requires data to span the entire solution domain. Solution need not begin or 
end at the table endpoints, however. 

Nonzero initial conditions must be specified for all variables in the IPINT 
array not internally computed as a function of flow parameters. 

This is done simply by using the commands VYY and VYYEND. As with other float- 
ing point data, the VYY Command card may also contain multipliers expressed as 
floating point numbers or integers pointing to locations in RARRAY. A negative 
sign (-) on an integer indicates multiplication by the reciprocal. Since no 
internal nondimensionalizing is performed on these arrays, the scale factors 
are usually used for this purpose. The VYYEND card must contain the variable 
number being initialized. For example: 

VYY -27 

210 * 300. T 

VYYEND 1 

Variable number 1 in the IPINT array is initialized to a velocity of 300 ft/sec, 
nondimensional ized by U located in RARRAY (27). 

oo 

In addition to variables in the IPINT array, other arrays may be initialized 
using the READ command. Starting address for storage in these arrays are stored 
in IZ array locations given in Appendix E. For example, initialization of nodal 
values of specific heat (see Appendix E), is accomplished by specifying variable 78 
on the READ command card: 

READ 5 300 78 0 

300 * 0.24 T 

The READ command above is interpreted as: READ 300 values from unit 5 and store 
them at the address specified in I Z( 78 ) plus 0. Similarly, data may be 
written from an array into a file for print, punch for later interrogation using 
the WRITE command. 

Solution Procedure (CMC Deck Section VII) 

C0M0C employs a flexible solution scheme which, through command card speci- 
fication, permits the user to call functional modules which perform sequences 
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of operations pertinent to a particular set of equations. The heart of the 
system lies mainly in the 5 link subroutines listed in Volume III. 

The 3DPNS equations utilize many of these, such as: Thermodynamics (LINK2-9), 
Turbulence Models (LINK5-6), Nonhomogeneous Laplacian Equation (LINK2-7), and 
Gradient Boundary Conditions (LINK1-7), to name a few. The order of solution 
for the 3DPNS equations proceeds by first calling DIMEN (LINK3-4) to compute 
dimensional izing parameters. This is followed by a call to GEOMFL (LINK1-3), 
to generate the finite element data. Nonlinear coefficients in the equations 
being integrated are evaluated at each integration step. The command LINKCALL 
is utilized to input the sequence of subroutine calls to be issued at each step, 
and the data are input in the form (link No., call list). For example, a data 
specification of: 

LINKCALL -1 

51, 29, 5 6, T 

would cause the program to sequentially call the subroutines (L INK5 ( 1 ) , 

LINK2(9), LINK5(6) prior to evaluating the derivatives at each integration 
step. By simply changing the linkcall list, therefore, one can vary the degree 
and type of nonlinearity in the solution. In the example, calling L INK5 ( 1 ) 
allows the geometry to vary according to the F21, F22 functions input under 
conmands VU3P0S and VU3VAL as previously described. 

TURBULENCE MODELING 

The 3DPNS algorithm is structured to accommodate a variety of turbulence 
closure models. The program presently addresses an algebraic mixing length model, 
and a more general turbulence kinetic energy (TKE) model, solving two 
differential equations written on TKE and dissipation functions (Vol. I, p. 11). 
While the models appear straightforward in concept, a variety of coefficients 
are required to support their application, thus complicating their usability. 

Table 2 defines the pertinent coefficients and suggested values for general 
3DPNS and boundary layer flows. The values are input as namelist data as 
indicated in the table and three principal options of use are available; 
mixing length theory, TKE, and a combination of both. Mixing length (MLT) is 
simplest to implement. Setting NE1E2 = 1 and ITKE = 0 will evaluate turbu- 
lent viscosity at each solution node until the solution has marched to the 
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Table 2 


Turbulent Viscosity Models 
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&NAME02 &NAME01 










station denoted by E1E2SW in (NAME02). The two equation (TKE-DISS) model is 
implemented by Getting the namelist values noted in the center column of Table 2. 

In most instances, initial values of TKE and DISS are not available. They 
can, however, be estimated using MLT. The input data for this option would 
appear as in column one of Table 2. Initialization is accomplished on turbulent 
assumption, and switches are preset in namelist NAME02 to initialize the TKE 
and DISS profiles using MLT. At each step below C4EDSW, TKE and DISS are eval- 
uated based upon the calculated mixing length. At station C4EDSW, ITKE is inter- 
nally set to one and (-NEQADD) equations are switched from parameter status to 
differential equation variables. This is internally accomplished by increasing 
the NEQADD counter by two. Utilization of this option requires that variables 
in the I PINT array be arranged such that the TKE & DISS variable numbers are at 
the initial NEQINT ( =NEQKNN & NEQADD) + 1 and + 2 locations. For 3DPNS, turbu- 
lent viscosities are evaluated in subroutine DFCFBL as coefficients for the 
viscous diffusion term in each of the momentum equations. 
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EXTERNAL JUNCTURE CORNER SAMPLE CASE DATA 


The Juncture region test case consists of turbulent flow parallel to a 
corner formed by the intersection of two 10 percent thick parabolic arcs with 
coincident leading edge. Figure 9 illustrates the flow field geometry. For 
parabolic Navier-Stokes, the solution is initiated at the origin, in Figure 9, 
over a generated grid on the secondary flow plane (z,y). The turbulent momen- 
tum equations with continuity constraint are subsequently integrated in the 
positive X direction (Fig. 9) and under the influence of the three-dimensional 
inviscid pressure field. Boundary conditions on the flow field include zero 
velocities along the walls and free-stream conditions at the outer plane. 

This section provides a detailed description of data requirements for the 
juncture region test case configuration, fully described in Volume I (section 
5). The data deck description parallels the decomposition noted in the previous 
section and concentrates on description of specific pertinent variables. Each 
subsection of the deck is listed with line numbers and footnote-like descrip- 
tions of each line which follow immediately. A continuous listing of the com- 
plete test case data deck is given in Appendix A. 
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II. 1. NAMELIST (Integer Scalar) Input 

1 FENAHE 

2 tNAMEOl 

3 NEQ ■ 9, NIZS - 250. NSNODE « 19. NSELEM « 2. ISIDE * 4. 

4 NVAR * 3. I PHI = 1. NODE * 193. LCOL » 100. KROU = 100. 

5 NDP = 10. NBAND « 23. LG * 32. NE1E2 * 1. NDERIV * 2. 

6 NDIM - 200. NHOUT « 2. NMBOUT * 50. NC » 8. NOUTS « 1. 

7 NOUTVC * 10. IBL * 1. NODES ■ 170. NPVSX « 2640. NHDL = 8. 

8 NEQAV2 * 1. NEQAV3 - 1. NPTS * 101. JPR =* 1. 

9 INOCRN - 1. IDDXST - 2. NTABPT = 14. NLINE = 140. 

10 NUALLS * 4. ICORN * 1. NPVSXT * 10, NNROW = 9, 

11 IULSEP - 19. NEQKNN = 5. NEQADD = -At NIHPLT = 1. 

12 KNTPAS = 200. NCNTIT * 0, NCNADD * 8. 

13 LPSUP * 6. LPPNCH ■ 2. 

14 SEND 

Line Command 


1 FENAME 

2 &NAME01 

3 NEQ 
NSNODE 
NSELEM 
ISIDE 

4 NVAR 

NODE 

LCOL 

KROW 

5 NDP 
NBAND 

LG 

NE1E2 


6 NMOUT 
NMBOUT 
NC 

7 NODES 
NPVSX 

8 JPR 


Command to initiate Namelist read 

Integer Namelist Data Array Name 

Number of dependent variable and parameter arrays 

Number of subdomain grid points specified 

Number of subdomains specified 

Maximum number of subdomain sides 

Number of variables (including coordinates) input 
for distribution over the refined grid 

Slightly larger than the number of nodes in the 
solution 

Larger than the number of nodes along abscissa 
(x 3 ) coordinate 

Larger than the number of nodes along the ordinate 
(x 2 ) coordinates 

One greater than NEQ 

Maximum bandwidth of Jacobian Matrix 

Number of columns in solution field 

0 + No Mixing length theory (MLT) used 

1 MLT used for diffusion coef. 

2 Delay using MLT until E1E2SW is set to 1 

(3,2) Print output in (Geometric form. Tabular form) 

Number of variables to be printed 

Number of characters in print fields 

Greater than maximum subdomain generated 
node density 

Number of pressures in P Versus Xj table 

1 indicates that geometric progression ratios 
are input to form a nonuniform grid 
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Line Command 

9 NLINE 

11 NEQKNN 

NEQADD 

12 KNTPAS 
NCNADD 

14 SEND 


Lines output per page control 

Number of ordinary differential equations solved = 
NEQKNN + NEOADD 

Number of ordinary differential equations 
initially not solved (Negative) 

Maximum number of integration steps between prints 

Number of integration steps prior to integrating 
secondary flow velocity equations 

NAMELIST End 


II. 2. NAMELIST (Real Scalar) Input 


15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 


INAME02 

UINF * 
DSTART 
TMULT * 
XMUINF 
XNUGEO 
CHI EPS 
ESCF = 
OSUSQ = 
AOMGEX 
RHSCAL 
CHITST 
HSINIT 

IEND 


100. Of TOFINF * 530. Of PINF = 2116. 8r REFL « I.Of 
« I.Of E1E2SU * 10000. Of COHPX « 2. Of COHPY = I.Of 
1 • 21 f TKEDGE « I.Of PIBAR * 2. Of PALPH * 20. Of 
- O.Of OSG = I.Of ALC = I.Of BLTH * I.Of EFhULT = I.Of 
» I.f C4EDSW - .001 f DELMLT = 1.E-4f PCFACT = I.f 
a 3.E-4F RHOIh « I.F PRDIS = 1.3f GUhULT * 1 . f 
I.Of RUEDSW * I.Of XLAh * O.Hf UCMULT * I.f OLDMLT = l.i 
1.0E-5f OHEGXP =* 1 . 5 f XDELTA = O.OIf EPSMIN = 1.0E-5f 
« 2. 0 f GAMFAC * 1.0E-20f GAhEXP = 9. Of BEXP * 4. Of 
« I.Of PPFACT * I.Of U2STRS = I.Of DELP * 101. Of 
= 10. Of CHISTP * 5. Of TO = O.OIIf SIMPLT ■ O.OIIf 
= 1 . OE-4 f TD * 0.589 f HMAX ■ I.Of TSADD * 0.05f 


Line 


Command 


15 

&NAME02 

Real Namelist Data Array Name 

16 

UINF 

Reference velocity 


TOFINF 

Reference temperature 


PINF 

Reference pressure level 


REFL 

Reference length 

17 

DSTART 

Initial step-size multiplier 


E1E2SW 

Station at which to reset NE1E2 


COMPX 

Geometric form print x 3 compression factor 


COMPY 

Geometric form print x 2 compression factor 

18 

TMULT 

Step size multiplier 

19 

XMUINF 

Reference viscosity 


ALC 

Characteristic Finite Element Size 
(minimum side length) 


EFMULT 

Maximum TKE initial level (ABS ValUe) 
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Command 


C4EDSW 

DELMLT 

CHIEPS 

RHOIM 

PRDIS 

ESCF 

RUEDSW 

XLAM 

VEMULT 

VLDMLT 

OSUSQ 

OMEGXP 

XDELTA 

EPSMIN 

AOMGEX 

GAMFAC 

GAMEXP 

BEXP 

RHSCAL 

PPFACT 

U2STRS 

DELP 

CHITST 

CHISTP 

TO 

HSINIT 

TD 

HMAX 

TSADD 

SEND 


Primary flow coordinate at which turbulence 
differential equations evaluation begins 

Turbulence equations source term scale factor 

ODE integration convergence factor 

Viscous wall damping factor 

Prandtl number for dissipation function 

Turbulence length scale factor 

Boundary layer measure 

Turbulence coefficient 

Secondary velocity convection term specifier 

Secondary velocity laminar diffusion 

Minimum TKE level 

Turbulence parameter 

Convergence criteria factor 

Minimum diffusion level factor 

Exponent on wall damping factor 

Diffusion factor 

Diffusion factor 

Diffusion factor 

Additional Reynolds stress turbulence model terms 

Implements the perturbation pressure in 
secondary flow velocity equations 

Implements Reynolds stress terms in secondary 
flow equations 

Print interval (% of total primary flow coord. Dim.) 
Maximum iteration count before step-size decrease 
Minimum iteration count before step size increase 
Initial primary flow direction coordinate 
Initial integration step size 
Primary flow coordinate total length 
Maximum step size (% of TD) 

Implements additional perturbation pressure terms 



1 1. 3 DYNAMIC ARRAY DIMENSIONING AND DEPENDENT VARIABLE SPECIFICATION 


29 FEDIMN 

30 IPINT -1 

31 15623789 0 Or 0 0 0 7*0» 10*11 1 T 


Line Command 

29 FEDIMN Dimension arrays to fit problem size 

30 IPINT Cards following specify dependent variable 

and parameter arrays 

Variable I primary flow velocity 

2 Secondary flow velocity (vertical) 

3 Secondary flow velocity (horizontal) 

5 Turbulence kinetic energy 

6 Dissipation function 

7 Perturbation pressure (p p ) 

8 Complementary pressure (p c ) 

9 Continuity equation potential (<h) 


Dependent variables and parameter arrays are 
stored in sequential arrays each of which is 
NODE long 


In NAME01 

NDP is the number of arrays for space allocation 

NEQ is the number of equations to be integrated 

NEQADD is the number of equations to begin inte- 
grating following initialization 
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III. FLOW-FIELD GEOMETRY, Nondimensionalizatiofivqnd Fiofte. Element Matricres. ; 


32 LINK4 9 

33 IARRAY 450 250 

34 PSIBD -1 

35 1 -2 T 

36 IARRAY 61 0 

37 LINK2 14 

38 NETA 

39 9 9 T 

40 NEPS 

41 9 9 T 

42 STYPE 

43 2*4 T 

44 SELCN 

45 9 10 3 1 12 13 19 2f 

46 1 3 14 15 19 16 17 18 T 

47 DEPVAR 289 290 1248 T TRANS • COORD* » NORMAL COORD * r U1 VEL ♦ 

48 0.0 0.869565 0.01918 5*0.0 0.0999 0.0999 

49 0.0 1.20 .869565 .01918 0.0 1.15 1.20 1.15 1.20 

50 0.0 0.0 .01918 5*0.0 0.0 .01918 

51 3*0.0 .0999 .0999 4*0.0 


52 

53 

0.0 

DONE 

0.0 . 

758 

5*0.0 0.0 .783 3*0.0 

.783 0.0 4*0.0 T 

54 

READ 

5 

63 

-26 


55 

11 

1 2f 

2 

12 11 T TURN DIAGONAL 

IN LOWER RIGHT CORNER. 

56 

READ 

5 

63 

-26 918 



57 181 171 182 1 171 172 182 T TURN DIAGONAL IN UPPER LEFT CORNER. 

58 READ 5 63 -26 972 

59 110 100 90 T ADD TRIABLE TO OUTSIDE CORNER. 

60 IARRAY 14 325 T ADD 1 TO ELEMENT COUNT. - NELEM * 325 


Line 

Command 


32 

LINK4 9 

Dynamically Dimensions Discretizer Arrays 

33 

IARRAY 

Reset location 450 in the scalar integer array 
to the value 250 

34 

PSIBD 

Set diagonals, each subdomain 

37 

LINK2 14 

Initiates discretizer 

38 

NETA 

Number of dimensions in the local coordinate 
system direction for each subdomain 

40 

NEPS 

Coordinate system direction for each subdomain 

Number of divisions in the e local coordinate 
system direction for each subdomain 

42 

STYPE 

Number of corner verticies on each subdomain 


33 



Line 


Command 


44 SELCN Grid point numbering for subdomain local 

coordinate definition 

1 . Order i s : 

Corner grid points counter-clockwise followed 
by side grid points counter-clockwise 

2. Corner grid points 1 and 2 define the local 
abscissa direction 

Corner grid points 1 and 4 define the local 
ordinate direction 

3. Nodes are generated in rows sweeping from the 
origin to the maximum ordinate 


45 


53 


55 


60 


DEPVAR Numbers on this card define the relative addresses 
for storage of coarse grid data following: 

289 indicates abscissa coordinate values 

290 indicates ordinate coordinate values 

1248 indicates primary flow velocity data 

DONE Signifies end of discretizer data and causes the 

grid refinement to occur and primary flow velocity 
values to be distributed over the refined grid 

READ Reads integer data from unit 5 (card reader) into 

array beginning at relative address 63 (generated 
finite element numbers). The data are displaced 
in the array 918 words from the beginning 

IARRAY Reset the generated finite element count to 325 


In NAME01 


NSELEM is the number of subdomains specified 

NSNODE is the number of specified grid points 

ISIDE is the maximum number of sides 

NODES dimensionalizes discretizer generated node- 
length arrays. Must be greater than - the 
maximum number of nodes to be generated in 
a subdomain 

NVAR is the number of arrays to be parametrically 
distributed 


34 


IV. TITLES AND HEADINGS 


Ol 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 
101 
102 


ue.at.Kxri 

WING / FUSELAGE JUNCTURE FLOW* 

DONE 

COMTITLE 

WING / FUSELAGE JUNCTURE FLOW* 

DONE 

DESCRIPT 332 T IOPAR PARAMETER TITLES FOR OUTPUT. 


REFERENCE 

ENGLISH-FT ENGLISH-IN M- 

-K-S C- 

-G-S 



1 FNRTH . ....... 

.FT 

T M 


.CM. . . 



un nr tty 


U A 

.M/S 

♦CM/S. 



DENSITY 

•LBM/FT3. • • • 

•N. A •••••••• 

KG/M3 ....... 

G/CC.. 


4 ♦ 

TEMPERATURE. * . 

• RANKINE • • * • 

•N. A ........ 

.KELVIN 

♦ N ♦ A. • 


• 444 

ENTHALPY 

. BTU/LBM. * * * 

♦N. A. •«**..» 

.KJ/KG...... 

. N . A . . 


4 4 4 4 

FROZ. SPEC. HEAT 

• BTU/LBM-R. • 


.KJ/KG-K.... 

•N.A. » 


4 4 4 4 

VISCOSITY 

• LBM/FT-S* • ♦ 


.NT-S/M2. ... 

♦POISE. . 

4 4#4 

LOCAL PRESSURE 

♦ PSF 


KNT/M2 

♦TORR. 

. . 

4 4 4 4 

LOCAL SOLUTION 

•MACH NO..** 

. DPDX1 

.ENERGY 

.INT. 

VAR. . . 

T^PFIX . i * » . » » • 


# (322 


i F 1 i i i i 



wunrnM u • q . . . . 

t t ..... . 
♦H31 ••♦.*•*♦ 






Xl/LREF ....«.• 

• DX1/LREF * * . 

.EPSILON... . 

.DX1M/LREF • • 

♦ REFL 

RE 

NO. 


DONE 

MPARA -1 

5*2 f 2 2 162 164 163f 3*2 164 163f 3*2 170 174f 
3*2 165 2 » 2 -175 3*2f 3*2 176 2f 3*2 177 178f 
2 2 169 168 167 f 3*2 108 2f 5*2f 5*2f 5*2 T 

IONUMB -1 

999 f 5*200 f 999 f 200 4*43f 200 27 200 27 27f 

200 10 200 10 10f 200 58 200 58 200f 

200 97 200 97 200f 200 30 200 30 200 f 

200 38 200 38 38f 999f 39 4*36f 300 154 100 135 122f 

398 4*11 186 f 200 4*11 139f 11 12 14 A5 47 T 

DESCRIPT 203 T TITLES FOR OUTPUT DEPENDENT VARIABLES • 


Ul/UREF 

U2/UREF 

U3/UREF 

NU/NUREF 

TKE/TKER 

BISS/D I SSREF 

PP / PSTAG 

U' 

U'V' 

U'U' 

V' 

DONE 

IOSAVE -1 

V'W' 

W' 




1248 2248 3248 1247 5248 f 6248 7248 3271 4271 5271 f 
6271 7271 8271 T 
IOMULT -1 

3*2 21 3*2 -222-22 -2f 13*1 T 


Line Command 

61 DESCRIPT 204 Std print titles 

65 COMTITLE Problem identifying titles 

08 DESCRIPT 332 Std output header labels 

£3 MPARA Scalar parameter print multipliers 

87 IONUMB Scalar parameter print locations 

93 DESCRIPT 203 Dependent variable print labels 

98 IOSAVE Dependent variable print locations* 

101 IOMULT Dependent variable print multipliers (RARRAY LOC.) 


*Note in line 99 that the dependent variables are identified by decoding each 
number specified (i.e. 1248 is stored in the first NODE locations beginning at 
the address in IZ( 248) ; Variable 2, 2248, is directly behind Variable 1 since 
it is stored in the second NODE locations beginning at the address stored in 
IZ ( 248) ; etc.). 


35 



V.l BOUNDARY CONDITIONS 


162 

IBORD -1 


163 

10*11 It 

17*110 20. 10*1-1 190. 18*1-10 171 T 

164 

KBNO 1 


165 

ADD 

DONE 

1 66 

19*110 1 

T FIX WAUL NODES. 

167 

KBNO 2 


168 

ADD 

DONE 

169 

19*110 1 

T FIX WALL NODES. 

170 

KBNO 3 


171 

ADD 

DONE 

172 

19*110 1 

T FIX WALL NODES. 

173 

KBNO 7 


174 

ADD 

DONE 

175 

19*110 1 

T FIX WALL ONLY FOR PP 

176 

KBNO 8 


177 

ADD 

DONE 

178 

19*110 10 

T 

179 

KBNO 9 


180 

ADD 

DONE 

181 

9*110 10f 

9*110 110 T FIX OUTSIDE BOUNDARY - BUT NOT CORNER. 

182 


WHICH IS NOW INTERIOR. 

Line 

Command 


162 

IBORD 

List of adjacent domain boundary nodes 

164 

KBNO 1 

Nodes where Ui is held constant (wall nodes) 

167 

KBNO 2 

Nodes where u 2 is held constant (wall nodes) 

170 

KBNO 3 

Nodes where u 3 is held constant (wall nodes) 

173 

KBNO 7 

Nodes where perturbation pressure held at intial state 

176 

KBNO 8 

Nodes where complementary pressure held at initial state 

179 

KBNO 9 

Nodes where equation potential (<(>) is held at zero 


Named-side cards are formatted [4(A8, 18, 4X)] 

Boundary conditions are applied to secondary flow plane. 


Positive direction is counter- 
clockwise (ie. [-RIGHT 2] fixes 
the right column and node numbers. 
Node numbers are stored beginning 
at the second row from the top and 
proceeding toward the bottom row). 



left 



right 


x 3 
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V. 2 INITIAL CONDITIONS 


183 

IARRAY 

47 

19 





184 

LINK3 

4 






185 

I ARRAY 

250 

0 

T I PHI IN GEOMFL 




186 

LINK1 

3 






187 

LINK2 

30 

T 

NODPCP 




188 

CNTPTS 

-1 






18? 

19*10 

T 






190 

CNTNDS 

-1 






191 

190*11 

1 

r 





192 

I ARRAY 

47 

19 





193 

LINK! 

2 






194 

LINK1 

11 

T 

NODPPR 




195 

LINK2 

10 

T 

TBLINP 




196 

LINK1 

11 

T 

NODPPR 




197 

IARRAY 

374 

190 




198 

RARRAY 

390 

1* 

0 




199 

LINK2 

3 






200 

RARRAY 

63 

1.0 




201 

LINK2 

15 






202 

LINK2 

3 






203 

LINKS 

6 






204 

LINK2 

15 






205 

LINK2 

3 






206 

LINKS 

6 






207 

RARRAY 

95 

0.9E-3 




208 

LINKS 

4 






209 

IARRAY 

107 

0 

T SET NE1E2 = 0 




210 

KBNO 

5 






211 

ADD 



DONE 




212 

19*110 

lr 

19*110 10 T FIX WALL 

AND 

FREESTREAM 

NODES. 

213 

KBNO 

6 






214 

ADD 



DONE 




215 

19*110 

If 

19*110 10 T FIX WALL 

AND 

FREESTREAM 

NODES* 


Line 

Command 


184 

LINK3 4 

Compute the nondimer.sionalizing parameters 

186 

LINK1 3 

Compute the finite element matrix coefficients 

187 

LINK2 30 

Compute the complementary pressure from Cp table 

188 

CNTPTS 

Reset the nodal subset string lengths 

190 

CNTNDS 

Reset the nodal subset strings 

193 

LINK1 2 

Generate finite element-node global connection table 

194 

LINK1 11 

Evaluate initial distributed complementary pressure (P c ) 

195 

LINK2 10 

Generate the primary flow velocities 

196 

LINK1 11 

Re-evaluate initial distributed complementary pressure (P c ) 

199 

LINK2 3 

Evaluate wall skin friction effects 

201 

LINK2 15 

Evaluate effective boundary layer integral parameters 

202 

LINK2 3 

Evaluate wall skin friction effects 

203 

LINK5 6 

Evaluate initial turbulence diffusion coefficients 
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V.I PRIMARY FLOW PARAMETER TABLES 


103 CNTPTS -1 

104 19 19 r 

105 CNTNDS -1 

106 10*1-10 100. 9*1-1 9f 10*110 100. 9*1-1 189 T 

107 I ARRAY 47 2 

108 LINK2 23 T CALL CPSTUP TO READ IN CP DATA* 

10? 0.0 0.02784 0.04731 0,07165 0.10210 0.14010 0.19000 

110 0.26870 0.41840 0.76770 1.51600 3.01300 6.00800 9.00200 

111 0*00999 

112 0.46940 0.46940 0.42770 0.38980 0.36030 0.33690 0.31740 

113 0.29980 0.28560 0.28160 0.28580 0.29330 0.33690 0.29710 

114 0.01914 

115 0.38290 0.38290 0.35650 0.32150 0.29170 0.26730 0.24680 

116 0,22830 0.21360 0.20940 0*21440 0.21980 0.24730 0.22400 

117 0.03126 

118 0.29870 0.29870 0*28500 0.25710 0.22970 0.20600 0.18560 

119 0.16700 0.15260 0*14910 0.15500 0.16010 0.17520 0.16460 

120 0.04629 

121 0.21490 0.21490 0.21020 0.19090 0.16830 0.14700 0.12800 

122 0.11040 0.09699 0.09488 0.10170 0.10680 0.11750 0.11180 

123 0.06409 

124 0.13190 0.13190 0.13350 0.12220 0.10560 0.08826 0.07195 

125 0.05642 0.04503 0.04494 0.05295 0.05821 0.06638 0.06353 

126 o.OR4«;fl 

127 0^05079 0.05079 0*05659 0*05191 ^."04182 0.02944 0.01691 

128 0*00491 -0.00370 -0.00108 0.00823 0.01375 0.02055 0.01944 

129 0.13310 

130 -0.10500 -0.10500 -0.09411 -0.08947 -0.08799 -0.08917 -0.09175 

131 -0.09444 -0.09397 -0.08425 -0.07182 -0.06549 -0.06010 -0.05904 

132 0.1609 

133 -0.1788 -0.1788 -0.1664 -0.1584 -0.1522 -0.148 -0*1451 0. 

134 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1*0 0. 

135 0.2224 

136 -0*3151 -0.3151 -0*3012 -0.2879 -0*2739 -0.2601 -0.2466 .2 

137 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 0. 

138 0.2908 

139 -0*4296 -0.4296 -0.4155 -0.3982 -0.3783 -0.3571 -0.3345 0 

140 -1*0 -1*0 -1.0 -1*0 -1.0 -1.0 -1.0 0* 

141 0.40250 

142 -0.54460 -0.54460 -0.53120 -0.51030 -0.48480 -0.45640 -0.42500 

143 -0.38670 -0.34020 -0.29330 -0.26800 -0.25890 -0.25390 -0.25000 

144 0.50 

145 -0.57430 -0.57430 -0.56140 -0.53960 -0.51270 -0.48250 -0.44880 

146 -0.40740 -0.35680 -0.30640 -0.28000 -0.27070 -0.26570 -0.26170 

147 0.5975 

148 -.5446 -.5446 -.5312 -.5103 -.4848 -.4567 -.425 

149 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 

150 0.71 

151 -.4296 -.4296 -.4155 -.3982 -.3783 -.357 -.3345 .7 

152 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 0. 

153 T 

154 VU2P0S -3 

155 0.01 0,4 1.0 T 

156 VU2VAL 136 

157 1,0 1.0 1.0 T 

158 VU3P0S -3 

159 0.01 0.4 1.0 T 

160 VU3VAL 136 

161 1.0 1.0 1.0 T 

Line Command 

103 CNTPTS Nodal subset string lengths 

105 CNTNDS Nodal subset strings 

108 LINK2 23 Reads pressure coefficient table (Cp) 

154 VU2P0S Primary flow coordinates 

156 VU2VAL X 2 coordinate variation at primary flow coordinates 

158 VU3P0S Primary flow coordinates 

160 VU3VAL X 3 coordinate variation at primary flow coordinates 
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in in 



VII 


SOLUTION PROCEDURE 


216 

LINKCALL -1 

217 

2 4 T 

218 

QKNINT 

219 

EXIT 

220 

CASE END 


CONTES DUMMY CALL* 


Line 

Command 


216 

LINKCALL 

List called at the end of each integration step 

218 

QKNNT 

Initialize implicit modified Newton Integration 

219 

EXIT 

Terminate run and return to MAIN for next data case 

220 

CASE END 

Last data case 
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APPENDIX A 

Juncture Region Data Deck 


1 

2 

3 

4 

5 

6 
7 

e 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


FENAME 

SNAME01 

NEQ * 9t NIZS * 250 f NSNODE « 19t NSELEH = 2r ISIDE « 4 ? 

NVAR * 3? IPHI « If NODE - 193f LCOL «= 100f KROU » 100r 
NDP * 10r NBAND * 23f LG * 32f NE1E2 * If NDERIV * 2f 
NDIM « 200 f NMOUT ■ 2f NMBOUT « 50f NC *= 8f NOUTS * If 
NOUTVC - 10f IBL * If NODES * 170f NPVSX » 2660f NNDL - 8f 
NEQAV2 * If NEQAV3 * If NPTS * 101 f JPR * If 
INOCRN * If IDDXST - 2f NTABPT = 14f NLINE - 140f 
NWALLS * 4f ICORN * If NPOSXT * 10f NNROW = 9f 
IWLSEP - 19f NEQKNN * 5f NEQADD * -4f NIhPLT * If 
KNTPAS * 200 f NCNTIT * Of NCNADD ■ 8f 
LPSUP - 6 * LPPNCH « 2f 

SEND 

SNAHE02 

UINF * 100. Of TOFINF « 530. Of PINF = 2116. 8f REFL * l.Of 
DSTART * l.Of E1E2SW * 10000. Of COMPX « 2. Of CONPY = l.Of 
TNULT « 1 .21 f TKEDGE « l.Of PIBAR = 2. Of PALPH » 20. Of 
XHUINF - O.Of OSG * l.Of ALC ■ l.Of BLTH «= l.Of EFHULT « l.Of 
XNWGEO * l.f C4EDSW = .001 f DELMLT » l.E-4f PCFACT = l.f 
CHIEPS * 3.E“4f RHQIM ** 1 . r PRDIS » 1.3f GUHULT » l.f 
ESCF « l.Of RUEDSW * l.Of XLAN * 0.11 f UCMULT « l.f ULDMLT * l.f 
OSUSG * 1.0E-5f OHEGXP » 1.5f XDELTA * O.Olf EPSMIN = 1.0E-5f 
AOHGEX « 2. Of GAHFAC * 1.0E-20f GAMEXP * 9. Of BEXP = 4. Of 
RHSCAL - l.Of PPFACT ■ l.Of U2STRS = l.Of DELP * 101 .Of 
CHITST * lO.Of CHISTP « 5.0f TO = O.Ollf SIMPLT «= O.Ollf 
HSINIT * 1 *0E-4 f TD « 0.589f HMAX - l.Of TSADD * 0.05f 

SEND 
FEDIMN 
IPINT -1 

156237890 Of 000 7*0f 10*11 1 T 

LINK4 9 

I ARRAY 450 250 

PSIBD -1 

1 -2 T 

IARRAY 61 0 

LINK2 14 

NETA 

9 9 T 

NEPS 

9 9 T 

STYPE 
2*4 T 
SELCN 

9 10 3 1 12 13 19 2f 

1 3 14 15 19 16 17 18 T 

DEPVAR 289 290 1248 T TRANS. COORD. f NORMAL COORD. f U1 VEL. 

0.0 0.869565 0.01918 5*0.0 0.0999 0.0999 

0.0 1.20 .869565 .01918 0.0 1.15 1.20 1.15 1.20 

0.0 0.0 .01918 5*0.0 0.0 .01918 

3*0.0 .0999 .0999 4*0.0 


0*0 

o 

o 

♦ 

758 5*0*0 

0.0 ♦ 783 

3*0.0 

♦783 0.0 

4*0.0 T 

DONE 







READ 

5 

63 -26 





11 

1 2f 

2 12 11 

T TURN 

DIAGONAL 

IN LOWER 

RIGHT CORNER. 

READ 

5 

63 -26 

918 




181 

171 

182 f 171 

172 182 

T TURN 

DIAGONAL 

IN UPPER LEFT CORNER. 

READ 


5 63 

-26 972 




110 

100 

90 T ADD 

TRIAGLE TO 

OUTSIDE 

CORNER. 



IARRAY 14 325 T ADD 1 TO ELEMENT COUNT. - NELEM - 325 
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*1 ue.t>UKiri 

62 WING / FUSELAGE JUNCTURE FLOW. 

63 

64 DONE 

65 COMTITLE 

66 WING / FUSELAGE JUNCTURE FLOW. 

67 DONE 

68 DESCRIPT 332 T IOPAR PARAMETER TITLES FOR OUTPUT. 


69 

REFERENCE 

ENGLISH-FT 

ENGLISH-IN 


M-K-S 

C-G-S 



70 

LENGTH 

.ft 

. , IN. » » ♦ 


. *M* * ....... 

♦ .CM. . . 



71 

VELOCITY 

.FT/S 

. *N. A. . • 



. .CM/S. 


♦ . * ♦ 

72 

DENSITY 

•LBM/FT3* . * 

. *N » A* . . 

. . . . 

t 

. KG/M3 

. G/CC.. 


• • 

73 

TEMPERATURE.... 

.RANKINE... 

. . N » A • • • 

♦ • ♦ • 

. .KELVIN.... 

. .N.A.. 


.... 

74 

ENTHALPY 

.BTU/LBM... 

♦ * N . A * . . 

. . * . 

. .KJ/KG 

. .N.A.. 


♦ 444 

75 

FROZ. SPEC. HEAT. 

• BTU/LBM-R. 

• *N. A. • * 

* • * • 

. .KJ/KG-K... 

. . N » A . * 

* * 

4 4 4 4 

76 

VISCOSITY 

• LBM/FT-S • . 

* • N . A. . . 

. ♦ . . 

. . NT-S/M2 . • • 

. .POISE.. 

4 4 4 4 

77 

LOCAL PRESSURE. 

. PSF 

. .PSI... 

* ♦ ♦ • 

. KNT/M2 

. .TORR. 

. . 

4 4 4 # 

78 

LOCAL SOLUTION. 

.MACH NO... 

. • DPDX1 . 

♦ • ♦ ♦ 

. .ENERGY.... 

. .INT. 

VAR ... 

79 

T2PFIX 

. H21 

. *G22... 


. fi?3 

• FI ... . 



80 

NWGEOM H 7 S ♦ . ♦ . ♦ 

,H31 




. .033 

• *G1 . . . 



81 

Xl/LREF 

• DX1/LREF. ♦ 

. .EPSILON... 

. .DX1M/LREF * 

. .REFL 

RE 

NO. 

82 

DONE 








83 

MPARA -1 








84 

5*2* 2 2 162 

164 163* 3*2 164 163* 3*2 

170 

174* 




85 

3*2 165 2* 2 

-175 3*2* 3*2 

176 2* 3*2 

177 

178* 




86 

2 2 169 168 167* 3*2 108 2* 

5*2* 5*2, 

5*2 T 




87 

IONUMB -1 








88 

999* 5*200* 

999, 200 4*43* 

200 27 200 

27 

27* 




89 

200 10 200 10 

10* 200 58 200 58 200* 






90 

200 97 200 97 

200, 200 30 200 30 200* 






91 

200 38 200 38 38* 999* 39 

4*36* 300 

154 

100 135 122* 




92 

398 4*11 186* 

200 4*11 139* 

11 12 14 

85 

47 T 




93 

DESCRIPT 203 

T TITLES FOR 

OUTPUT DEPENDENT 

VARIABLES. 




94 

Ul/UREF 

U2/UREF 

U3/UREF 


NU/NUREF 

TKE/TKER 



95 

DISS/DISSREF 

FP / PSTAG 

U 7 


U 7 V 7 

U 7 U 7 



96 

V 7 

V 7 W 7 

W 7 







97 DONE 

98 IOSAUE -1 

99 1248 2248 3248 1247 5248* 6248 7248 3271 4271 5271* 


100 

6271 7271 

8271 T 






101 

IOMULT -1 







102 

3*2 21 3*2 -222 

-2 2 -2* 13*1 T 




103 

CNTPTS -1 







104 

19 19 T 







105 

CNTNDS 

1 






106 

10*1-10 100* 9*1-1 

9* 10*110 

100. 9*1-1 189 T 



107 

I ARRAY 47 

2 






108 

LINK2 

23 T 

CALL CPSTUP 

TO READ 

IN CP DATA. 



109 

0.0 

0.02784 

0.04731 

0.07165 

0.10210 

0.14010 

0.19000 

110 

0.26870 

0.41840 

0.76770 

1.51600 

3.01300 

6.00800 

9.00200 

111 

0.00999 







112 

0.46940 

0.46940 

0.42770 

0.38980 

0.36030 

0.33690 

0.31740 

113 

0.29980 

0.28560 

0.28160 

0.28580 

0.29330 

0.33690 

0.29710 

114 

0.01914 







115 

0.38290 

0.38290 

0.35650 

0.32150 

0.29170 

0.26730 

0.24680 

116 

0.22830 

0.21360 

0.20940 

0.21440 

0.21980 

0.24730 

0.22400 

117 

0.03126 







118 

0.29870 

0.29870 

0.28500 

0.25710 

0.22970 

0.20600 

0.18560 

119 

0. 16700 

0.15260 

0.14910 

0.15500 

0.16010 

0.17520 

0.16460 

120 

0.04628 







121 

0.21490 

0.21490 

0.21020 

0.19090 

0.16830 

0.14700 

0.12800 

1 22 

0.11040 

0.09699 

0.09488 

0.10170 

0.10680 

0.11750 

0.11180 

123 

0.06409 







124 

0. 13190 

0.13190 

0.13350 

0.12220 

0.10560 

0.08826 

0.07195 

125 

0.05642 

0.04503 

0.04494 

0.05295 

0.05821 

0.06638 

0.06353 


126 0.08458 
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127 

0.05079 

0.05079 

0.05659 

0.05191 

- .“04182 

0.0294. 

0.01691 

128 

0.00491 

-0.00370 

-0.00108 

0.00823 

0.01375 

“0.02055 

0.01944 

129 

130 

0.13310 

-0.10500 

-0.10500 

-0.09411 

-0.08947 

-0.08799 

-0.08917 

-0.09175 

131 

-0.09444 

-0.09397 

-0.08425 

-0.07182 

-0.06549 

-0.04010 

-0.05904 


132 0 . 160 ? 


133 

-0.1788 

-0.1788 

-0.1664 

-0.1584 

-0.1522 

-0.148 

-0.1451 

0. 

134 

135 

-1.0 

0.2224 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

0. 

136 

-0.3151 

-0.3151 

-0.3012 

-0.2879 

-0.2739 

-0.2601 

-0.2466 

.2 

137 

138 

-1.0 

0.2908 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

0. 

139 

-0.4296 

-0.4296 

-0.4155 

-0.3982 

-0.3783 

-0.3571 

-0.3345 

0 

140 

141 

-1.0 

0.40250 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

0. 


142 -0*54460 -0.54460 -0.53120 -0.51030 -0.48480 -0.45640 -0.42500 

143 -0.38670 -0.34020 -0.29330 -0.26800 -0.25890 -0.25390 -0.25000 

144 0.50 

145 -0.57430 -0.57430 -0.56140 -0.53960 -0.51270 -0.48250 -0.44880 

146 -0.40740 -0.35680 -0.30640 -0.28000 -0.27070 -0.26570 -0.26170 

147 0.5975 

148 -.5446 -.5446 -.5312 -.5103 -.4848 -.4567 -.425 .5 

149 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 .5 

150 0.71 

151 -.4296 -.4296 -.4155 -.3982 -.3783 -.357 -.3345 .7 

152 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 -1.0 0. 

153 T 

154 VU2P0S -3 

155 0.01 0.4 1.0 T 

156 VU2VAL 136 

157 1.0 1.0 1.0 T 

158 VU3P0S -3 

159 0.01 0.4 1.0 T 

160 UU3UAL 136 

161 1.0 1.0 1.0 T 

162 IBORD -1 

163 10*11 If 17*110 20 f 10*1-1 190? 18*1-10 171 T 

164 KBNO 1 

165 ADD DONE 

166 19*110 1 T FIX WALL NODES. 

167 KBNO 2 

168 ADD DONE 

169 19*110 1 T FIX WALL NODES. 

170 KBNO 3 

171 ADD DONE 

172 19*110 1 T FIX WALL NODES. 

173 KBNO 7 

174 ADD DONE 

175 19*110 1 T FIX WALL ONLY FOR PP 

176 KBNO 8 

177 ADD DONE 

178 19*110 10 T 

179 KBNO 9 

180 ADD DONE 

181 9*110 10? 9*110 110 T FIX OUTSIDE BOUNDARY - BUT NOT CORNER? 

182 WHICH IS NOW INTERIOR. 

183 I ARRAY 47 19 

184 LINK3 4 

185 v I ARRAY 250 0 T IPHI IN GEOMFL 

186 LINK1 3 

187 LINK2 30 T NODPCP 

188 CNTPTS -1 

189 19*10 T 

190 CNTNDS -1 

191 190*11 1 T 

192 IARRAY 47 19 
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193 

LINK1 

2 






194 

LINK1 

11 

T NODPPR 





195 

LINK2 

10 

T TBLINP 





196 

LINK1 

11 

T NODPPR 





197 

IARRAY 

374 

190 





198 

RARRAY 

390 

1.0 





199 

LINK2 

3 






200 

RARRAY 

63 1,0 





201 

LINK2 

15 






202 

LINK2 

3 






203 

LINK5 

6 






204 

LINK2 

15 






205 

LINK2 

3 






206 

LINK5 

6 






207 

RARRAY 

95 

0.9E-3 





208 

LINKS 

4 






209 

IARRAY 

107 

0 T SET NE1E2 » 

0 




210 

KBNO 

5 






211 

ADD 


DONE 





212 

19*110 

If 

19*110 10 T FIX 

WALL 

AND 

FREESTREAM 

NODES 

213 

KBNO 

6 






214 

ADD 


DONE 





215 

19*110 

If 

19*110 10 T FIX 

WALL 

AND 

FREESTREAM 

NODES 

216 

LINKCALL 

-1 






217 

2 4 

T 

CONTES DUMMY CALL. 





218 

QKNINT 







219 

EXIT 







220 

CASE END 








APPENDIX B 


The CMC discretizer possesses the capability of refining coarsely dis- 
cretized flow domains having boundaries of general shape. By subdividing the 
flow domain into a few shapes having quadratic variation along the boundaries, 
virtually any closed domain can be accurately modeled and discretized to con- 
form to solution requirements. The method is graphically presented in Fig. 

B.l for a single element two-dimensional airfoil in an infinite domain. In 
Fig. B.la the flowfield has been divided into twelve subdomains consisting 
of eight quadrilateral and four triangular shapes. The grid points represent 
the data required to specify the subdomain boundaries. Note that all boun- 
daries are of degree 2 or less. The side grid points serve to describe the 
curved surface shapes and their location relative to adjacent verticies deter- 
mines the spacial distribution of the generated grid. Figure B.lb illustrates 
a 648 element grid generated from the subdomain description in Fig. B.la. 

Note that the refined grid around the airfoil surface and toward the leading 
and trailing edges is due to the off center placement of the subdomain side 
grid points. Variations on refinement type and number of elements is readily 
accomplished from a single set of subdomain data by simply moving side grid 
points and respecifying the number of elements for each subdomain. 

The grid refinement algorithm proceeds by looping over the subdomain and 
generating grid point and finite element correction data. For each of the 
quadrilateral and triangular subdomain shapes, a coordinate transformation 
exists which maps a local natural coordinate system onto the physical plane. 

The general form of the required transformation for all space is: 

x i = Qj x j j = 1, n (A.l) 

where n is the total number of specified grid points on each subdomain boundary 
and summing over repeated indices is observed. The scalar consists of the 
grid point physical values of the coordinates at the ith generated node. The 
quantity Qj is called the shape function and contains functions of e and n 
which satisfy the transformation illustrated in Figure B.2 for e and n of the 
{xj}. Equation (B.l) is valid over one-, and two-, and three-dimensional _ 
Cartesion space and shape factors may be derived for a variety of geometric 
shapes and polynomial degrees. The Q. for biquadratic functional shapes over a 
general quadrilateral for instance is J 


Q 3 


'4(1 - 
4(1 + 
4(1 + 

J 4(1 - 
4(1 - 
4(1 + 
4(1 - 
, 4(1 - 


e)(l - n)(-e - n - 1) 
e) (1 - e)(e - n - 1) 
e) (1 + n)(e + n - 1) 
e)(l + n)(-e + n - 1) 
e 2 )(i - n) 
e)(l - n 2 ) 
e 2 ) (1 + n) 
e)(l - n 2 ) 


IT 


> 


( 6 . 2 ) 
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where the e, n axis bisects the quadrilateral and the equations .are .ordered 
according to the gr'i d point numbering noted in Figure A. 2. Substitution of 
equation B.2 into B.l for a specific set of {x j } and evaluating the equation 
over the limits of c and n (-1 to 1) yields a biquadratic approximation 
to x over the subdomain. Accuracy of the values of x-j are dependent upon 
the ability of the shape functions to approximate the physical geometry. A 
curvature which is exactly biquadratic, for instance, will be interpolated 
exactly using equation B.2. Note in Figure B2b that the side located grid 
points lie at exactly midside, by definition. The side nodes in Figure B.2a, 
however, need not be at exactly mid-side since it is not required in the 
definition of {Q>. Movement of the specification of the side nodes in the 
physical plane, therefore, allows for a smoothly varying distribution of 
generated data over the subdomain as illustrated in Figure B.l. The coor- 
dinate transformation is of the serendipity family and can yield nonunique 
distributions if excessive skewing of the mid-side node is attempted; there- 
fore, care must be taken when applying the method as in Figure B.l, where 
regions of highly refined grid are necessary along the airfoil surface. 

Each of the subdomains of which the entire solution domain is composed 
may be treated independently if the data generated along subdomain boundaries 
is identical for each subdomain. This requires that intersubdomain boun- 
daries share the same grid points and that (e;n) be independent along the sub- 
domain interfaces. To illustrate that this is indeed the case, consider the 
two coincident subdomains of Fig. R.3. 

Along the common boundary (3,4,7), and (4,1,8) 2 


rij = C = 1.0 

£2 s c = -1.0 

and evaluating equation (2) with these constraints yields: 

Q * LTsCe + e 2 ), %(-£ + e 2 ), (1 - e 2 ) (B- 3) 

Q = LJs(-n + n 2 ), *s(n + n 2 ), (1 - n 2 ) (B. 4) 

And since e,® -n 2 along the boundary, equation (A. 3) becomes 

Q = lh (-n + n 2 ), %(n + n 2 ), (1 - n 2 ) ( B, 5) 
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Fig. B .3 Two adjacent subdomains . 


which is identical to equation (B.4). Other boundary combinations and the tri- 
angular shaped subdomain can be shown to yield similar results. Three- 
dimensional subdomains have surface interfaces and similar equations can be 
derived as above where the equivalent of equation (B.5) contains functional 
expressions fora planar surface. This demonstrated generality, therefore, 
completely eliminates any requirements for specifying the subdomains in a 
particular order and numbering may begin at any vertex. 

It is impossible to encompass the entire spectrum of geometric shape and 
refinement requirements.within a few subdomain types. The completeness of 
a data refinement scheme, therefore, requires that addition to its geometric 
flexibility be open ended. Utilizing the subdomain method, special subdomain 
functions may be developed and inserted directly into the system. An example 
of this is a telescoping quadrilateral having quadratic interpolation. The 
purpose of this special function is to provide a means for increasing or de- 
creasing the number of generated triangular elements when progressing from 
one end of a subdomain to another. Fig. B.4 illustrates the effect for a 
quadrilateral expanding from three elements per side to six elements per side. 
This subdomain has proven very useful for flows having small areas of high 
velocity gradient where a very fine mesh is required relative to the overall 
flowfield mesh size, and windowing-in is not practical. Fig. B.5 illustrates 
such a case where vortex generation at a wing tip is to be computed. Very 
high velocity gradients occur in the immediate area of the tip and decay 
rapidly. Utilizing the telescoping subdomain, the grid size is reduced from 
six elements per side in the tip region to two elements per side where an 
infinite boundary condition was applied and 64 elements were eliminated from 
the computational domain. 

The interpolation functions described above are used to distribute the 
grid point data over each subdomain. Combining the subdomains to form the 
complete discretization, nowever, poses the problem of having duplicate sets 
of data along the interconnecting boundaries. Removal of the duplication is 
accomplished through the construction of a subdomain connection table from 
the specified data and interrogating it from an algorithm which operates over 
the subdomain boundaries. Table B.l illustrates the connection table for the 
first four subdomains of Fig. B.J. The storage requirements for the table 
are 2l(no. of subdomain sides). 
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TABLE B*1 Subdomain Connection Table for Joukowski Airfoil. (First 4 Subdomains) 


Subdomain 

Side 

Subdomain 

Side 

1 

1 

10 

3 

1 

2 

0 

0 

1 

3 

2 

1 

1 

4 

0 

0 

2 

1 

1 

3 

2 

2 

0 

0 

2 

3 

3 

1 

2 

4 

0 

0 

3 

1 

2 

3 

3 

2 

0 

0 

3 

3 

4 

1 

3 

4 

0 

0 

4 

1 

3 

3 

4 

2 

0 

0 

4 

3 

5 

1 

4 

4 

0 

0 


Note that the first two columns need not be stored since they are sequential. 

For the solution domain of Fig.B .1, therefore, only eighty words of storage 
are required to administrate boundary duplicatiorr elimination. The algorithm 
proceeds subdomain by subdomain, first generating a sequence of local 
(dummy) node numbers, then interrogating the connection table to determine 
if any sides are connected to any other subdomain. Upon finding one which 
is of lower number it retrieves the local boundary node numbers for that sub- 
domain and substitutes them for the coincident set into the proper locations. 
Simultaneously, the duplicate generated variable set at the boundary is elim- 
inated. The finite element connection table is subsequently formed for all 
generated elements in the subdomain and sequential global node numbering is 
substituted in a simple loop at the end. The connection table is also use- 
ful for locating external boundaries since they appear as zeros in the table. 
Note in the above scheme that only the subdomain data were manipulated to 
sequence the generated data thus maintaining higher efficiency of the algorithm. 
Also, storage requirements are minimized by requiring operation on not more 
than two subdomains simultaneously, hence the data for a large problem could 
be generated interactively on a minicomputer or small time sharing system 
utilizing disk or tape intermediate storage. An especially useful feature of 
such a method when combined with a graphics package and a video (CRT) terminal 
is the ability to interactively modify and debug data decks. The above scheme is 
very efficient and generated data for 24 variables over 10 nodes/CPU sec on 
an IBM 360-65- 
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Command Name Input 


A C CN T ROL CARD WITH THE FOLLOWING PARAMETERS IS READ I N - 


PARAMETER FORMAT CARD COLS. 


DESCRIPTION 


VI *3 1-3 CONTROL VARIABLE. 

NMOL fpcc AF r El? 8 NMULI1) = NX 

MMULI2) = NPCD 
NMULI3) = NREPET 
NMULI4) = NR T AP£ 


TO NMUL(N) 

IF KOIJMP = 1 !N ^r A *1- 01. THEN THE Emtjrp lupjj CARD IS PRINTED 
IMMEDIATELY AFT=R RFjNjr, P f a P FOLLOWED BY THE DATA THAT IS BEING 
STORED ALONG W! T H the Ca T a*< f N t R y POSITION IN T HE IZ ARRAY 




ARSVAL Nl N2 '13 

SE T PZ(!Z(N2J) = ABSIRZI IZIN3) ) « I = 1, Nl 

( BLANK > 

tr.r^b J T] SC^M ANOTHER CARO. 

cr^ic 

PRINT T H T COMOC TITLE PAGE TWO TIMES. 

cov t i t l: 

RFAO A TITLE CARD WHICH WILL APPEAR ON COMO". 

OLSC.PI^ NX 

call DSC pp>T AND PROCESS ACCORDING TO NX. 


NX 


BLANK - s E AD AND WRI^E INFORMATION CARDS. 


203 - READ TITLES FOR DEP. VAR. OUTPUT HEADINGS. 

204 - READ DESCRIPTIVE TITLE FOR HEADING 


A T BEGINNING OF OUTPUT HEADER 


332 - 5 EAD TITLES FOR PARAMETERS PPINTED IN 


THE 0UTO ! J T HEADER. 

FNC 

u'tTURN CONTROL TO MAIN PROGRAM, RESET ARRAYS 
AND c FT URN TO BDINPT. 

r XI~ 

call c xit. 

2DPF 

SE T 1 3 L = ) AND I PHI = 1. 

3 r b p 

S-t 1 3 L = 1. 

3Hp\)$ 

SET IBL = l AND N 3 DP NS = 1. 

2 n NS 

S- T IBL = 0. 

F c 0!MN 

CALL D IMENS I ON AL I Z AT I ON ROUTINE FLDIMN. 

FFMA^E 

CALL f = naME TO SET DEFAULT SCALARS AND THEN 
CALL f.MFLST T 0 READ IN NAMEOl AND NAME02 

namelis t ? . 
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DESCRIPTION 


CGMMANO N AM? 


IAFPAY Nit N2, Vi, N4, -TC. 

SET IARPAY(Ml) = N2t IARPAY(N3I * N4, ETC. 

ICCNH CALL IC.PND TO PRINT REAL AND INTEGER SCALARS. 

INDEXSET NI N2 rj3 N4 NS N6 M7 

DO FDR I = l, NX 
K = I Z ( I Z(N2)«-N3*I) - 1 
F Z C I Z ( NA ) «-N5 I = Ri ( IZIN6 ) *N7M ) 

INPl'T Nl SET INPUT UN I J TP Nl . 

INTEGER ALLOWS NEW VALUES TO READ INTO A SEQUENCE 

CF LOCATIONS IN THE BORDER, IPLACt AND NLOC 
VICTOR S. 

KBiN'"’ NX F N ' r c R FIXED MODES FDR DEP. VARIABLE NX. 

CALL. GET BNP 

KBN'l MX l ENTER FIXED NODES AND/OR BOUNDARY 

CONDI “ !CNS FOR D c P. VARIABLE NX. 

CALL GE T RCD 

LINK l NX CALL LINK l (NX) 

LINK? NX CALL L INK 2 ( N X ) 

LTNK3 NX CALL L INK3 ( NX , DUMMYl , DUMMY2 > 

L I N< A NX CALL LINKAUX.K) 

LINKS NX CALL L I NK 5 ( N X ) 

MATSU‘S Nl, N2» N3, NA, NS 

CALL N A T SUM ( PZ(IZ(N2)J, RZUZ(N3)>, RAPE AY ( NA) , 
?Z(IZ(N5>), Nl ) 

N A M - L 1ST CALL KM e lST T 0 READ IN NAMEOl AND NAME02 

NAYFLISTS. 

PDUMP M N2 M3 CALL °DUMP ( IZ(IZ(N2)), IZ(IZ(N3)), NL ) 

PLUS M M2 M3 . . . ET . 

R ARP AY INI) = ?ARFAY(N2) * R ARR A Y ( N 3 ) ♦ ...«■ 

kAROAY Nl N2 . . . Cl . . . C.N . . . NN 

RARRAY(Nl) * PAPRAYIN2 ) * . . . Cl * . . . R ARRAY 

P c AD Ml M2 N 3 NA 

R EAT ( Nl ) P Z ( 1 1 ( M3 )+NA ) , I = 1 , M2 

FETPIEV- Nl M2 CALL NBNOPY (1, IZ(N2>, IZ(IYY) ) WITH NP = Nl. 

R c ST~.p Ml N2 CALL F.ES T 0F ( Nl, I Z ( N2 ) , IZ(IYY) ) 

SAV6TAPC Ml N^ SAVr OUTPUT PN UNIT NL, REWIND AFTER 

MOD ( KCUNT, N2 I = 0. 

S F 7 V A L M N? NT NA NS 

CALL SETVAL ( RZUZIN2)), RZUZIN3II, 

PAPFAY(NA), R ARRAY ( N5 ) , Nl ) 

RQRT Nl N2 ... PARPAY(Nl) = S QRT ( R ARRA Y ( N 2 ) 

PAP.RAY(N3) = SQRT (RARRAY (NA) ) 
c T C * 

SQRT Nl N 2 N 3 RZ(IZ(N2t = S QRT ( R Z ( I Z ( N3 ) ) ) , I = 1, Nl 

VECMUL Nl N2 N3 NA 

DC FOP I = 0 TO NX-1 

FZIN2MI = P Z ( N 3 <■ I ) * RZ(NAH) 

VYY4ND NX DENOTES 'END °< r INPUT FOR DfcP. VAR. NX. 

WRITE Nl M2 N3 NA 

f WRITE ( Nl ) D Z( IZ(N3)+N4) , I = 1, N2 


INPUT Ml 
I N T E G E F 


KBN D 


L I NK l 
LINK? 
LINK 3 
LINK A 
LINKS 


NAM- L 1ST 

PDUMP M N2 N3 
PLUS M N2 M3 , 

k ARR AY N l N2 . 
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T H: FOLLOW I NG COMMANDS AP P UTILIZED TO LOCATE STORAGE ADDRESSES 
FOR IN T EGER ARRAY CATA SP = C IF IC AT TCN. THE CARDS ARE DIRECTLY 
FOLLOWED BY INTEGER DATA CARDS IN FREE FORMAT. THE READ SCAN IS 
TERMINATED BY a 'T' or blank card. 


IZUPLACEUn = LOCATION IN TH p IZ ARRAY AT WHICH TO BEGIN 
STORING INTEGER c N T R. I E S • 

IA°RAY(NlOC(K) ) = NUMBER 0= ENTRIES STORED STARTING AT 

I Z (I P L AC r ( K ) ) . 

IF NX ,NF. -I, ENTER L I T ER AL CATA. * 

I F NX . E Q . -I, EN T = o NUMERICAL CATA DIRECTLY. 


■ 

command 

NAME 

I 0 L AC c ( K ) 

NL 0 C ( K ) 

DESCRIPTION. 

3 

THICK 

63 

Q3 

ELEMENT THICKNESS VECTOR. 

7 

I PI NT 

6 

31 

SOLUTION SEQUENCE VECTOR. 

12 

LINKCALL 

121 

125 

LINK NOS. TO BE CALLED AT END OF QKNUI N. 

13 

SPEC IES 

31 

12 l 

VAPIA3LF NDS • FOR SECIES TO BE 3JN. 

14 

IOMijl t 

122 

67 

OUTPUT VARIABLE MULTIPLIER FROM RAPPAY. 

15 

T 0 s A V p 

124 

60 

VARIABLE LIST TO BE DISPLAYED AT JJTPUT. 

16 

CNTPIC 

12 7 

47 

CONTOUR NODES TO BE USED IN 
CONTES, DFCF3l, TRBTHK , WLFLXS, ETC. 

17 

CN^NDS 

123 

128 

NO. OF NODES IN EACH CONTOUR LINE. 

18 

IB33C 

33 

131 

C’OUNT-R-CLOCKWI SE LIST OF BOUNDARY. 

21 

i n nu ^ 13 

131 

142 

LIST OF ENTRIES IN RARRAY TO BE 
ntSPLAYFD AT START OF EACH 3UTPJT. 

22 

MpARA 

135 

67 

LIST CF MULTIPLIERS IN RAPPAY USED 
TO MULTIPLY IONUMB ENTRIES. 

26 

NX 

74 

67 

DUMMY 

27 

NY 

758 

67 

dummy 

28 

3 Lt NT $ 

26 

67 

PE AD IN ELEMENT NODE CONNECTIONS. 


* K L TT 

t 1 $u Rc> ar 

’INF •GET3N0*. 

N 0 ^ = - 

'NLGC(K) OR NVOCIKI = 67 IMPLIES 

P C E SF T LENGTH IS NOT CHANGED. 
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THE FALLOWING COMMANDS ARE utilized to locate storage aodrfsses 
FOR PEAL ARRAY DATA S p ECI F ICATICN. THE NAME CARDS ARE DIRECTLY 
FOLLOWED BY REAL DATA CARDS IN FRFfc FORMAT . T H£ READ SCAM IS 
’■ c R M I N A T r D BY A *T* OR BLANK CARD. 


IZ(NPLACFIK) ) = LOCATION IN T HE IZ ARRAY AT WHICH TO BEGIN 
STORING REAI ENTRIES. 

IARRAYINVuCK) ) = N'JM R C R Cj F r N TR I E S ST°RED STARTING A T 

!Z(NPLACr.(K) ) . 

T H f S E OCTROI CARDS CAN CONTAIN A GROUP OF MULTIPLIERS FOR 
thf entered oata. 

E . G . 

V Y Y 3 -IOJ.O -2 7 

MUL'IPLIES c A CH INRUT VECTOR NUMBER BY 

-100.0 * RARRAYI3) / KARP. AY ( 27 J 


K 

^ n *4 
NAMf 

N°l AC =( K ) 

NV’^C ( K ) 

n E$CR IPTIOM. 

t> 

VU3PCS 

65 

178 

X STATION for variable grid CHANGE. 

4 

V J?VAL 

65 

67 

SCALE F AC r OR FOR VARIABLE GRID CHANGE. 

5 

V T H I C K 

7 ) 

69 

V A LU c - OF ELEMENT THICKNESSES. 

default = i.o z alc 

6 

V*H n 

84 

67 

density node points. 

DEFAULT = PHOINF 

7 


19 

59 

TABLE LOOK-JP TEMPERATURES. 
default = tlfinf 

8 

VCP T AB 

1 3 

67 

table look-up specific heats. 

PE R At.!L T = CPJINF 

0 

VX ICOP 

8 9 

16 

XI -COORD I NATES AT NODc POINTS. 

10 

vxncr* 

9 0 

16 

X2-CP ORDINATES AT NODE POINTS. 

1 1 

VH 

79 

67 

ENTHALPY DISTRIBUTION AT NODE POINTS. 
HE FAULT = 1.0 

16 

VP ?: ss 

91 

6 7 

PR F SSUP.E VALUES AT NODE POINTS. 
DEFAULT = PINF 

IT 

VSCHMIO T 

114 

67 

SCHMIDT NO. DIST. AT NODc POINTS. 
DEFAULT = SCT 

18 ! 

V Y V 

3 2 

67 

Of. PENDENT VAR. DIST. A T NODE POINTS. 

19 

V Tr MP 

8 5 

6 7 

T =MO"P ATURE DIST. AT NODE POINTS. 

default = tofinf 

?p 

vtk 

8 8 

67 

THICKNESS n F ELEMENTS IN T HICK VECTOR. 
DEFAULT = 1.0 / ALC 


vsuthld 

l 

6 7 

$" r LD VR , SILDTR, STLDCRt STLDLX, STLCON 
ENTRIES FOR SUTHc? LANDS LAW. 

D C F. . 1 1 63E-4 » 4 C ;4.0, 204.0, 1 .5, 0.0 


54 








APPENDIX C ( Cont'd.) 


D 

COMMAND 

NPLACE(K) 

NVOC(K) 

DESCRIPTION 

24 

VPPANDTL 

134 

67 

PFANOTL NO. 0 1 ST. AT NODE POINTS. 
DEFAULT = pp 

23 

VX3ST 

139 

16 L 

DOWNSTREAM STATIONS AT WHICH PRESSJRE 
IS DEFINED. 

26 

v^vsx 

140 

67 

DOWNSTREAM PRESSURES AT VX3S^. 

27 

y/cpc: t L9N 

136 

67 

TU»B'JLENT VISCOSITY AT NODE POINTS. 
DEFAULT = XMJINF 

29 

r &r - a y 

0 

67 

F A&R AY ( NX ) = AMULT, WHERE 

AMULT = COMBINATION DF RcMAINING ENTRIES. 

31 

V W ALLS T A 

6 7 

32 

DOWNSTREAM PQS. AT WHICH TO 
INJECT TRANSVERSE VELOCITY. 

32 

VNAl LVAL 

6 3 

31 

VALUE OF INJECTFD TRANSVERSE VELOCITY. 
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Scalar Variable Descriptions 


VALUES IN PARENTHESIS ARE DEFAULT CONDITIONS. 


ALPHANUMERIC I ARRAY ENTRIES 


I ARRAY 

NAME 

ENTRY 

96 

I BASE 

261 

IBC 

99 

IBL 

344 

IBLAS 

207 

IBOT 

327 

IDDXST 

127 

IDIFRT 

162 

1EPSET 

347 

INCLFT 

348 

INCRGT 

158 

INITCN 

152 

INITKE 

295 

IPLOTV 

42 

INPUT 

28 

IPASS 

105 

IPTSPL 

132 

IPWRIT 

294 

IRAT 

100 

IREND 

259 

ISIDE 

322 

ISTART 

215 

ISUPRS 


186 

ITDA 

187 

ITtlB 

205 

I TIMER 

97 

I TKE 

133 

I TOP 

196 

ITWALL 

312 

IUONLY 

122 

IWRIT 

92 

IZSIZE 

291 

JPR 


169 KCDC 
61 KDUMP 
167 KNTPAS 


DESCRIPTION 


(200) RASE NO. FOR 12 ENTRIES. 

NUMBER OF BOUNDARY CONDITION TYPES. 

1 * BOUNDARY LAYER PROGRAM 

EPS EXPONENT FOR U2» U3 CONVERGENCE CRITERION. 

1ST ELEMENT AT UHICH TO PRINT DEBUG INFORMATION. 

< 1 > 

NO. OF TIMES TO PRINT INTER. OUTPUT IN LINK1 » GETPPR. 
JNCPPR AND PRSGRD. 

1 = FORCE COMPUTATION OF DISSIPATION FUNCTION. 

OUTPUT VECTOR INCREMENT FOR SETVAL AND MATSUM. 

1ST INPUT VECTOR INCREMENT FOR MATSUM. 

INPUT VECTOR INCREMENT FOR SETVAL. 

2ND INPUT VECTOR INCREMENT FOR MATSUM. 

INITIALIZER IN CONTES. 

1 = TKE. DISSIPATION ARE ALREADY INITIALIZED. 

DATA SET NUMBER ON UHICH TO STORE DATA FOR PLOTTING. 

(5) INPUT LOGICAL UNIT NUMBER. 

NO. OF CALLS TO DERVBL . 

(0) * USE LUDWIEG - TILLMAN FORMULA FOR TAU WALL. 

1 = USE PATANKER AND SPALDING'S FORMULA FOR TAU WALL. 
DEBUG CODE IN LINK3 AND STRF FOR INTERMEDIATE OUTPUT. 

SET = TO 1 WHEN N3DPNS IS FIRST TURNED ON. 

END POSITION IN 'IZ' ARRAY. 

NUMBER OF SIDES / SUPER ELEMENT. 

STARTUP CODE IN IMPLCT. 

1 = SUPPRESS PRINTOUT OF? 

A. STARTUP IN BDINPT . 

B. NODE MAP. 

C. OUTPUT STATIONS. 

UNIT NO. ON UHICH TO STORE INTEGRAL PARAMETER DATA 
FOR PLOTTING. 

UNIT NO. ON UHICH TO STORE 'PLOTS' DATA FOR PLOTTING. 

SET = NO. OF TIMES TO CALL TIMETK SUBROUTINE. 

NOTES TIMETK NEEDS TO BE DEFINED FOR YOUR INSTALLATION. 

0 = DO NOT INTEGRATE TKE - DISS. EQUATIONS. 

1 = USE TKE - DISS. TO COMPUTE TURBULENT VISCOSITY. 

2ND ELEMENT AT WHICH TO PRINT DEBUG INFORMATION. 

1 = USE DUDY FOR TAU WALL. 

(2) NO. OF STEPS UNTIL CONVERGENCE ON U1 ONLY. 

DEBUG PRINT FLAG IN CONTES. DFCFBL » GETFSL AND TAUU. 
MAXIMUM DIMENSION OF IZ VECTOR. 

1 - -MIDSIDE NODES FOR ETA DIRECTION ARE 
PROGRESSION RATIOS. 

1 * RESET NLINE TO 50 AND DUMP KODE TO 2. 

PRINT INPUT CARDS AND DATA GENERATED IN BDINPT. 

(99) MAXIMUM NO. OF INTEGRATION STEPS BETWEEN PRINTS. 
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6 KODG - PRINT GEOMETRY OUTPUT IF *NE* 0* 

7 K0D5 - PRINT INTERMEDIATE DERVBL OUTPUT K0D5 TIMES# 

26 KOUNT - RUNNING COUNT OF OUTPUT. < LIMITED BY LPRINT « ) 

113 KPLVAR - <10> NO* OF VARIABLES TO BE PLOTTED OR PUNCHED* 

86 KPNT - 1 = PRINT STATION ( SET DURING EXECUTION*) 

8 KPRINT - PRESENT VALUE OF PRINT COUNTER* 

52 KROU - (100) NO* OF ROWS IN DISCRETIZATION* 

151 KWFLXS - ( LCOL ) NO* OF SLICES AT WHICH TO COMPUTE TAU WALL* 

50 LCOL - (20) NO* OF COLUMNS IN DISCRETIZATION. 

47 LG - NO. OF COLS* IN SOLUTION FIELD. 

IF • NE * 0 ON INPUT * THEN CNTPTS AND CNTNDS ARE 
TO BE READ IN* 

330 LMDJAC - (1) COMPUTE JACOBIAN EACH ITERATION* 

214 LMLT - (LG) NO* OF CONTOURS FOR WHICH TO COMPUTE MIX. LENGTH 
TURBULENT VISCOSITY. 

179 LOC - INTERVAL NO. FOUND IN LOOK SUBROUTINE. 

172 LOWD - (2) USE LAMINAR VISC* BELOW LOWD AND MLT FROM LOWD ON* 

376 LPHI - +/- 1 TO EVALUATE PHI OR PP RHS IN PPRES* 

353 LPLOT - FLAG TO INITIATE PRINT OF DATA AT PLOT STATION. 

(SEE DPSISQ AND REOUTP) 

354 LPLTPR - FLAG TO ACTUATE PRINT OF DATA AT PLOT STATION (REOUTP) 
352 LPPNCH - NO* OF PASSES BEFORE LPUNIT TAKES AFFECT FOR K* EPS. 

34 LPRINT - (100) LIMIT ON OUTPUT COUNT. 

351 LPSUF* - NO. OF DATA STATIONS BEFORE LPLTPR TAKES AFFECT. 

350 LPUNIT - UNIT NO. ON WHICH TO STORE PPr Ul* U2? U3p Kf EPS DATA* 
212 MLTRHS - (1) NUMBER OF RIGHT HAND SIDES TO SOLVE FOR IN STRF • 

377 MONE 1 IN PPRES TO ALTERNATE SIGN OF LPHI. 

23 NB - (4) NO* OF CHAR. IN EACH WORD OF OUTPUT VAR* TITLE. 

305 NBAND - (2*LC0L+3> MAXIMUM BANDWIDTH OF JACOBIAN MATRIX. 

170 NBC - MAX. NO. OF BOUNDARY COND* FOR ANY ONE DEP* VAR. 

131 NBORD - NO. OF NODES AROUND BORDER OF DISCRETIZATION* 

22 NC - (8) NO. OF CHARACTERS IN OUTPUT FORMAT. 

125 NCALLS - (10) NO. OF ROUTINES TO CALL AT END OF INTEGRATION STEP 
108 NCNADD - BEGIN U2 * U3 INTEGRATION AFTER NCNADD INTEGRATION STEPS 

173 NCOMOC - NO. OF CARDS READ IN FOR COMOC TITLE PAGE. 

210 NCOMPG - (35) NO. OF ENTRIES TO PRINT FROM PRGDUM COMMON BLOCK* 

174 NCOMTD - NO* OF CARDS READ IN FOR TITLE INFORMATION. 

307 NCONV - NON-CONVERGENCE FLAG SET IN IMPSLV* 

13 NCOORD - FLAG FOR GENERATING AXI-SYMMETRIC DATA IN GEOMFL . 

59 NCPTAB - (1) NO* OF ENTRIES IN SPECIFIC HEAT TABLE* 

1 ND - INITIALIZATION PARAMETER 1WDFCFNS. 

317 NDBGPT - NO. OF ENTRIES TO PRINT IN VECTORS FOR DEBUGGING. 

124 NDERIV - (2) = CALL DERVBL • 

194 NDP - (10) SPACE ALLOCATION IN IIPINT VECTOR. 

14 NELEM - NUMBER OF ELEMENTS IN SOLUTION. 

304 NELPAS - INTEGRATION STEP COUNTER. 

31 NEQ - (5) MAXIMUM NO. OF VARIABLES TO BE SOLVED. 

43 NEQADD - NO. OF DIFF. EONS. INITIALLY NOT SOLVED. 

- E.G. -2 = DELAY INT* TKE AND DISS. 

UNTIL C4EDSW IS SATISFIED. 

286 NEQAV2 - WHEN NCNADD IS SATISFIEDt NEQADD = NEQAV2 + NEQAV3 . 

287 NEQAV3 - WHEN NCNADD IS SATISFIED* NEQADD = NEQAV2 + NEQAV3 ♦ 

58 NEQKNN - (1) NO. OF DEP. VAR* TO BE INTEG* IN GKNUIN* 

37 NEWPRT - (5) NO* OF SCALARS TO PRINT ACROSS A PAGE (MAX. « 8). 
107 NE1E2 - 0 = DO NOT USE MIXING LENGTH THEORY FOR DIFF. COEF . 

1 = USE MLT FOR SOLUTION OF DIFF* COEF. 

2 * DELAY USING MLT UNTIL E1E2SW IS SATISFIED. 

46 NF - (4) NO. OF y NB * BYTE WORDS IN TITLE FOR EACH DEP. VAR. 

54 NHHALF - PASS AT WHICH MAXIMUM ITERATIONS TOOK PLACE* 

53 NDB - CURRENT PASS (TOR PRINT POINT). 

68 NI - STARTING LOC. IN DEP* VAR ♦ MATRIX FOR THIS VARIABLE. 

301 NIMPLT - (1) IMPLICIT INTEGRATION SCHEME BEING USED* 



APPENDIX D (Cont'd.) 


I ARRAY NAME 
ENTRY 

316 NIT 


302 NITER 
94 NIZS 
323 NLAST 
88 NLINE 
191 NM 


306 NMBJAC 
60 NMBOUT 
206 NMDL 


190 NMOUT 

193 NM2 
16 NNODE 
340 NNROW 
55 NODE 
233 NODES 

19 NOE 
325 NOUEDG 
142 NOUTPR 

40 NOUTS 
38 NOUTVC 
30 NP 
157 NPASS 
336 NPDBUG 

198 NPGRDT 

199 NPGRDV 

20 NPRNT 

153 NPUNCH 
161 NPVSX 
371 NPVSXT 
313 NR 
328 NRJACB 

67 NS 

154 NSD 
252 NSELEM 
146 NSFDBE 

27 NSKIP 
64 NSM 
251 NSNODE 
189 NSTD 

155 NS2 
394 NTABPT 
140 NTCNTS 
203 NTCRDM 

62 NTITL 

197 NTPRNT 

177 NU2P0S 

178 NU3P0S 

260 NVAR 

90 NYY 

91 NZZ 
192 N2M 
166 N3DPNS 


DESCRIPTION 


KOBE FOR DEBUG PRINT OF: 

1* JACOBIAN 

2* L/U DECOMPOSITION 

ITERATION COUNTER ♦ 

(200) NO* OF IZ ENTRY POINTS THAT CAN BE STORED* 

KOBE FOR LAST STATION WHEN DX ♦ LT ♦ H. 

(60) LINE COUNT OUTPUT CONTROL • 

(3) TYPE OF ELEMENTS IN SOLUTION. 

2 = LINE (ONE-DIMENSIONAL) ♦ 

3 = TRIANGLE (TWO-DIMENSIONAL). 

<NBAND*NNODE) NO. OF ENTRIES IN JACOBIAN. 

(30) NO. OF VARIABLES TO BE PRINTED. 

(8) ALLOW EXTRA STORAGE IN IZ(71> AND IZ(72-76) 

LENGTH OF IZ(71) = MAXIMUM (NODE, NODE*NMDL) 

LENGTH OF IZ(72-76> = MAXIMUM ( NODE r (NODE/2 ) *NMDL ) 

(3) = PRINT OUTPUT IN GEOMETRY FORM. 

2 * PRINT OUTPUT IN NODE NO. SEQUENCE. 

NM**2. USED FOR STORING FULL MATRICES. 

NUMBER OF NODES IN SOLUTION. 

NO. OF COLS. (ROWS) ON ONE SIDE OF JUNCTURE CORNER. 

(100) VARIABLE DIMENSIONING PARAMETER IN FEDIMN. 

(NODE) VECTOR LENGTH FOR REFINE GRID GENERATOR. 

NO. OF EQUATIONS BEING SOLVED FOR DEP . VAR. 'NP'. 

FLAG TO STORE UEDGE FROM CP INFORMATION INTO U1 VECTOR. 
(100) NO. OF SCALARS TO PRINT IN OUTPUT. 

(10) NO. OF OUTPUT VECTORS TO PROCESS AT ONE TIME. 

(5) NO. OF OUTPUT SCALARS ACROSS PAGE AT PRINT STATION. 
DEP. VARIABLE BEING SOLVED AT THIS TIME. 

NO. OF PASSES THRU DERVBL ♦ 

NO. OF ITERATIONS FOR DEBUG IN DERVBL AT DEBUG POINT. 

(4) STARTUP COUNTER USED IN PRSGRD. 

(4) STARTUP COUNTER USED IN PRSGRD. 

(132) NO. OF PRINT POSITIONS ON A LINE OF OUTPUT. 

SET = 7 IF ELEMENTS AND NODES ARE TO BE PUNCHED IN DIMEN. 
(2) NO. OF PRESSURES IN P VS X TABLE. 

NO. OF ENTRIES IN EACH TABLE. (USED WITH NTABPT). 

NO. OF PRINTS FOR IMPLICIT INTEGRATION DEBUG. 

(1) COMPUTE JACOBIAN EACH ITERATION. 

GENERAL DUMMY PARAMETER. 

INITIALIZATION CODE IN DFCFBL . 

NO. OF SUPER ELEMENTS IN GRID GENERATOR. 

RESET CONDITION FLAG IN ' FINDBE ' . 

(NODE) NO. OF BOUNDARY LOC . / DEP. VAR. 

(10) STOP PROGRAM IF ANY OUTPUT EXP. IS .GT. NSM. 

NO. OF SUPER NODES IN GRID GENERATOR. 

(8) NO. OF STANDARD MATRICES TO BE STORED. 

INITIALIZATION CODE IN DFCFBL. 

NO. OF TABLES TO BE READ IN CPSTUP. 

STARTUP PARAMETER IN CONTES. 

STARTUP PARAMETER IN XYCRDM ♦ 

(10) NO. OF TITLE CARDS TO BE READ IN AND PRINTED AT THE 
BEGINNING OF EACH OUTPUT SET. 

99999 = DO NOT PRINT INTEGRAL PARAMETERS IN TRBTHK * 

(20) MAX. NO. OF ENTRIES FOR VAR. GEOMETRY DEFINITION. 
(CROSS PLANE) 

(20) MAX. NO. OF ENTRIES FOR VAR. GEOMETRY DEFINITION. 
(TRANSVERSE PLANE) 

NUMBER OF VARIABLES TO BE DISTRIBUTED OVER REFINED GRID. 
(4) SETS (NODE) / DEP. VAR. LOCATION IN IYY VECTOR. 

(4) SETS (NODE) / DEP. VAR. LOCATION IN IZZ VECTOR. 

(NM*2) USED FOR STORING SYMMETRIC MATRICES. 

SET * 1 WHEN JV .LE. ( NEQKNN+NEQADD ) IN IMPLCT. 
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RARRAY NAME 
ENTRY 


156 

AINF 

5 

AJ 

3 

ALC 

393 

AM8 

275 

AOMGEX 

86 

AVD 

137 

AVDP 

279 

BEXP 

209 

BLTH 

176 

CBTOKJ 

311 

CC1 

312 

CC2 

313 

CC3 

314 

CC4 

365 

CD 

356 

CE 

303 

CHIEFS 

302 

CHISTP 

322 

CHITST 

304 

CIMPTH 

364 

CK 

211 

CF0V2 

83 

COMPX 


84 

COMPY 

124 

CON 

70 

CONRHO 

62 

CONV 

207 

CONVRG 

309 

CPHI1 

310 

CPHI2 

158 

CPA 

159 

CPH 

160 

CPINF 

30 

CPOINF 

153 

CVCP 

148 

CVH 

151 

CVP 

152 

CVRHO 

150 

CVT 

149 

cvu 

357 

CW 

182 

C1D0RF 

183 

C1K0RE 

183 

C2D0RF 

184 

C2K0RE 

315 

C2C4 

143 

C4EDSW 

318 

DCHECK 

336 

DELCHK 

338 

DELMLT 

13 

DELP 

205 

DELSTR 

337 

DELTST 

103 

DEPLT 

165 

DRTODK 

319 

DSTART 

390 

DYNPRS 


DESCRIPTION 


- REFERENCE SPEED OF SOUND » 

- <778.28) JOULES CONSTANT. 

- < MIN * SIDE) CHARACTERISTIC ELEMENT SIZE. 

- TERM 8 MULTIPLIER IN PPRES FOR PP EQUATION. 

- (2.0) EXPONENT ON WALL DAMPING FACTOR 'OMEGA ' 

- (25.3) DAMPING FACTOR IN DFCFBL . 

AVDP = AVD * SQRT < RHO/RHOWAL ) * ANULOC / ENUT 

- < 4 * 0 > GAMMA = 0.01 * DELTAY * < Y/DELT A ) **BEXP 

- BOUNDARY LAYER THICKNESS* DELTA. 

- (4.184) SPECIFIC HEAT BRITISH TO MKS ♦ 

- COEF. FOR REYNOLD STRESSES COMPUTED FROM CPHI1 AND CPH2 

IN DIMEN. 

- SAME AS CC1. 

- SAME AS CC1. 

- SAME AS CC1. 

- (0.09) TKE - DISS. COEF. 

- <1.0) DIFF. COEF. MULTIPLIER FOR TKE. 

- (1.0E-4) CONVERGENCE FACTOR FOR IMPLICIT INTEGRATION. 

- < 4 . 0 ) MAX. NO. OF ITERATIONS FOR INCREASING STEP SIZE. 

- <10.0> MAX. NO. OF ITER. BEFORE DECREASING STEP SIZE. 

- (0.5) RELAXATION FACTOR FOR IMPLICIT INTEGRATION. 

- <1.0) TKE - DISS. COEF. 

- SKIN FRICTION 

- COMPRESSION FACTOR FOR OUTPUT COL. VECTOR 

INDICATES PERCENT OF X3 AXIS TO BE USED TO SHORTEN 
SPACING INTERVALS. 

- COMPRESSION FACTOR FOR OUTPUT ROW VECTOR. 

SAME AS COMPX, BUT FOR X2 AXIS. 

- <0.435) KARMANN'S CONSTANT USED IN MLT IN DFCFBL. 

- IF .GT. 0.0, SET ALL RHO = CONRHO. 

- <1.0>- OUTPUT SCALE FACTOR = 1.0 / REFL. 

- REAL NO. EQUIVALENT OF NCONV (NO. OF ITERATIONS). 

- (2.8) EMPIRICAL CONSTANT FOR COMPUTING CC1 - CC4 

- (0.45) EMPIRICAL CONSTANT FOR COMPUTING CC1 - CC4 

- (.24*32.174) SPECIFIC HEAT OF AIR. 

- (3.445*32.174) SPECIFIC HEAT OF HYDROGEN. 

- SPECIFIC HEAT COMPUTED IN CPINIT. 

- (0.24*31.174) REFERENCE SPEDIFIC HEAT. 

- (4186.0) SPEC. HEAT CONVERSION FACTOR. 

- (1.0) ENTHALPY CONVERSION FACTOR. 

- (4.725E-4) PRESSURE CONVERSION FACTOR. 

- (16.02) DENSITY CONVERSION FACTOR. 

- (1.0) TEMPERATURE CONVERSION FACTOR. 

- (0.3048) VELOCITY CONVERSION FACTOR. 

- (1.0) SCALE FACTOR FOR DISSIPATION DIFF. COEF. 

- (1.44) DISS. FCT • PRODUCTION TERM COEFFICIENT. 

- (1.0) TKE PRODUCTION TERM COEFFICIENT. 

- (1.92) DISS. FCT. DISSIPATION TERM COEFFICIENT. 

- (1.0) TKE DISSIPATION TERM COEFFICIENT. 

- ( CC2*CC4 ) COEF. FOR DIFFUSION COEFFICIENTS. 

- (30000.0) TKE - DISS. STARTUP POSITION IN DFCFBL. 

- ( DSTART *SSINIT)**3) 

- ( DELTST*DELMLT ) GEN. AND DISS. TERMS FOR TKE* EPS 

ARE SET TO 0.0 WHEN, DELUSQ .LT. DELCHK. 

- (1.0E-5) SCALE FACTOR FOR DELCHK. 

“(2.0) PERCENT INTERVAL FOR PRINTOUT* 

- DISPLACEMENT THICKNESS. 

- FIRST DERVBL PASS FOR LARGEST GRAD*U**2 ON ELEMENT. 

- (101.0) PERCENT OF TD TO BE USED FOR PLOTTING STATIONS. 

- (5. 0/9.0) DEGREES RANKINE TO DEGREES KELVIN 

- (10.0) SCALE FACTOR FOR DCHECK • 

- SET .GT. 0.0 FOR CONSTANT PRESSURE FIELD. 
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RARRAY NAME 
ENTRY 


175 EBTOKJ 
358 EFMULT 

90 EKNINF 
204 ENER 
135 ENERGY 
108 ENMULT 

14 EPS 

89 EPSINF 
274 EPSMIN 

95 EPTEST 
68 EP4MD 
368 ESCF 
171 EULER 
145 E1E2SW 

1 FACT 
80 FACTH 

59 FACTMU 
79 FACTP 

163 FTTOCM 
162 FTTOIN 

164 FTTOMT 
189 FI 
329 FI 0 

327 F2 
31 G 

277 GANEXP 
276 GAMFAC 

60 GAMJ1AF 
263 GSCALE 
353 GUMULT 
142 G1 
381 G10 

328 G2 

187 G22 

188 G23 

140 G32 

141 G33 

15 H 

16 HMAX 
45 HS 

7 HSINIT 
12 HT 
186 H21 

386 H21L 
139 H31 

387 H31L 
272 OhEGXP 

2 ONE 
396 OSG 
375 0SH1SQ 


DESCRIPTION 


(2.3244) ENTHALPY BRITISH TO MKS. 

(-0.01) ABS(EFMULT) BECOMES LARGEST LEVEL 
OF TKE INITIALIZED. 

(UINF**2> TKE NON-D FACTOR. 

DIMENSIONAL ENERGY FOR U1 VELOCITY . 

NON-DIM. ENERGY FOR U1 VELOCITY. 

DIMENSIONAL I ZING FACTOR FOR ENERGY. 

(XMUINF*RE*UINF*UINF/(G*ALC**(NM-1 > > 

CONVERGENCE FACTOR FOR INTEGRATION STEP. 

(UINF**3/ALC) NON-D FACTOR FOR DISSIPATION. 

(1.0E-5) SCALE FACTOR ON XhUINF FOR MINIMUM 
LEVEL OF DIFFUSION. 

* EPMULT « ZERO TEST FOR DISSIPATION USED IN DERVBL . 

- (1.0) MULTIPLIER FOR XMDOT IN PRSGRD. 

- (3.0) SCALE FACTOR IN TKE AND DISS. LENGTH. 

- <PSTAG/(RH0INF*UINF**2> EULER NUMBER. 

(30000.0) STATION AT WHICH TO SET NE1E2 = 3 - NE1E2. 

<ALC) NON-DIMEN ♦ FACTOR. 

(1.0 / (CP0INF*T0FINF) ) 

(RHOINF*UINF*ALC) 

(1.0/ FACTMU) 

- (30.48) FEET TO CENTIMETERS. 

- (12.0) FEET TO INCHES. 

- (0.3048) FEET TO METERS. 

- Y-COORDINATE OF FI CURVE. 

LAST VALUE OF FI CURVE. 

Y-COORDINATE OF F2 CURVE. 

(1.0) GRAVITATION CONSTANT. 

(9.0) GAMMA = 1.0/(1.0 + GAMFAC* ( Y/DELTA ) **GAMEXP ) 
(1.0E-20) SEE GAMEXP DEFINITION. 

(1.4V FACTOR USED IN GAS LAWS. 

(0.1) SCALE FACTOR FOR U1 MULT. OF GRADIENT ( PH I ) IN PPRES. 
KODE TO ADD GRADIENT ( PHI ) TO RHS OF U2> U3 EQNS . 
Z-COORDINATE OF G1 CURVE. 

LAST VALUE OF G1 . 

Z-COORDINATE OF G1 CURVE. 

VARIABLE GEOMETRY FACTOR. 

VARIABLE GEOMETRY SCALE FACTOR. 

- VARIABLE FACTOR. 

- VARIABLE FACTOR. 

- CURRENT TRIAL STEP SIZE. 

- (2.0) MAX. PERCENT OF TD TO USE AS STEP SIZE. 

- CURRENT STEP SIZE. 

- (1.0E-5) START INTEGRATION STEP SIZE AT THIS VALUE. 

OUTPUT VAR. FOR TIME STEP = HS * FACT / REFL 

(1.0) GRID GROWTH SCALE FACTOR. 

LAST VALUE OF H21. ( USED IN XYCRDM ) 

(1.0) GRID GROWTH SCALE FACTOR. 

LAST VALUE OF H31 . ( USED IN XYCRDM ) 

(1.5) TKEEXP « 2.0 * ( 2.0 - OMEGXP ) 

(1.0) PROGRAM CONSTANT. 

ADD STRESSES TO PP EQUATION IN PPRES. 

- (1.0 / H21**2> 
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RARRAY NAME 
ENTRY 

136 OSHAX 
271 OSUSQ 

362 0512 
361 0S6 

363 0S60 

339 PCFACT 
174 PDF TOC 
170 PDFTOK 

36 PEDDIM 
39 PEDGE 
268 PHICOD 
82 PIBAR 
9 PINF 
180 PMSKGS 

340 PPFACT 
99 PPRCON 

100 PPRIME 
67 PR 
185 PR D I S 

166 PSFTOA 
169 PSFTOI 
168 PSFTON 

167 VSFTOT 
351 PSTAG 

20 PTIM 
345 PUMULT 
179 RADCON 

21 RE 
43 REFL 
47 REFLRE 

285 RELAX 

289 RHOIM 
10 RHOINF 

157 RHOUIN 
360 RHSCAL 
199 RNULOC 

290 ROMULT 
119 ROUALC 

116 RR 

32 RTC0N1 
56 RTC0N5 

117 RT0HM1 


DESCRIPTION 


HEIGHT OF SOLUTION DOMAIN AT INITIAL STATION* 

< 1 ♦ O/TKEINF ) MINIMUM TKE LEVEL ALLOWED. 

1 • O/FACTORI AL < NM+1 ) 

1 * 0/FACT0RIAL<NM) 

< 2 . O/FACTORI AL < NM+2 > ) 

ADD PC TERM TO U2, U3 EQUATIONS. 

(0.01602) P0UNDS/FT**3 TO GRAMS/CM**3 
(16.02) P0UNDS/FT**3 TO KG/M**3 
DIMENSIONAL PRESSURE « PEDGE * PSTAG. 

NON-DIM. PRESSURE AT PRESENT STATION. 

(1.0/H) GRADIENT PHI MULTIPLIER AT END OF PPRES. 

UTAU COEF. FOR COMP. U1 IN JNCINP. 

FREESTREAM PRESSURE. DEF. = 1ST VALUE IN P VS X TABLE 
(1.0/2. 2) POUNDS / KG. 

USE PP TERM IN U2, U3 EQUATIONS. 

( RH0INF*UINF**2/ALC ) 

PRSSURE GRADIENT COMPUTED IN PRESSURE ROUTINE. 

(1.0) PRANDTL NUhBER ♦ 

(1.3) DISSIPATION PRANDTL NUMBER. 

( 4 ♦ 725E-4 ) P0UNDS/FT**2 TO PSIA 
(6.924E-3) P0UNDS/FT**2 TO P0UNDS/IN**2 
(47.88) P0UNDS/FT**2 TO NEWT0NS/M**2 
(0.3591) P0UNDS/FT**2 TO TORR. 

(PINF+0.5*RH0INF*UINF**2) 

PRINT TIME PARAMETER IN 'IMPLCT'. 

WHEN = 1, USE EDDY VISCOSITY FOR U2r U3 DIFF. TERM. 
(57.2957759) CONVERSION FACTOR RADIANS TO DEGREES. 

(RHOINF*UINF*ALC/XMUINF) REYNOLD'S NO. 

(1.0) REFERENCE LENGTH. 

(RHOINF*UINF*REFL/XMUINF> REYNOLD'S NO. BASED ON REFL. 
(0.6) RELAXATION FACTOR FOR GRID GROWTH IN LOOKAV. 

(1.0) USE WALL DAMPING FOR TKE* EPS EONS. IN DERVBL ♦ 
(PINF*XMA/<RUNIV*TSINF> FREESTREAM DENSITY. 

RHOINF * UINF 

ADDIT. TERMS FOR REYNOLD STRESSES IN RNLDST • 

1.0 - USE LOCAL VISCOSITY FOR VAN DRIEST DAMPING FACTOR. 
(PINF*XMA/<RHOINF*RUNIV*TOFINF) CONV. FACTOR IN DRHOBL. 
RHOINF * UINF * ALC**2 

CPH / CPA 

2.0 * G * AJ 

UINF**2 / <RTC0N1*CP0INF*T0FINF> 

RR $ ( TOH/TOA - 1.0 ) 
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RARRAY NAME 
ENTRY 

280 RUEDSU 

138 RUESO 
28 RUNI V 
129 SCT 
210 SHPFAC 
301 SIMPLT 
190 SLOPE 
50 SSINIT 
75 STLCON 

73 STLDCR 

74 STLDEX 
72 STLDTR 
71 STLDVR 

321 STPMLT 


374 TAREA 
372 TBAR 
35 TD 

22 TF 
206 THETA 

4 THK 

23 TIME 
48 TIMESV 
42 TKEDGE 
90 TKEINF 

300 TMNTS 
305 TMULT 

24 TO 

146 TOA 

58 TOFINF 

147 TOH 

40 TRATIO 
284 TSADD 
347 TSCALE 
155 TSINF 
111 TWELVE 
26 TWOPI 
394 T2FIX 
398 T2PFIX 
397 T3FIX 
371 UBAR 


DESCRIPTION 


- (10.0) IF ( Y/DELTA .GT. RUEDSU ) 

EPSDIM = GAMREF*UEDGE*DSTAR*GAMMA 

- RHOEDG * UEDG * UEDG 

- ( 1545 ♦ 33*32*174) UNIVERSAL GAS CONSTANT « 

- (1.0) CONSTANT SCHMIDT NUMBER* 

- SHAPE FACTOR* 

(TO) STATION AT WHICH TO START IMPLICIT INT* 

SLOPE OF VARIABLES COMPUTED IN LOOK* 

- (HSINIT / FACT) 

- WHEN *GT • 0*0* HAVE CONSTANT VISCOSITY* 

- (204.0) REF. CON* TEMP. IN SUTHERLAND. 

- (1*5) EXPONENT USED IN SUTHERLAND. 

- (492.0) REF. TEMP. USED IN SUTHERLAND. 

- (1.163E-5) VISCOSITY USED IN SUTHERLAND* 

- (1.1) STEP SIZE MULTIPLIER WHEN 

ALL VARIABLES CONVERGE ON 1ST ITERATION. 

- TOTAL COMPUTATIONAL AREA. 

- MASS WEIGHTED AVERAGE TEMPERATURE. 

- (1.0) TOTAL SOLUTION TIME < DISTANCE ) FROM TO. 

- (TF=TO+TD) FINAL STATION* 

- MOMENTUM THICKNESS. 

- (1.0) DEF. NON-DIM. THICKNESS OF ELEMENTS. 

- CURRENT STATION* 

- SAVED TIME LOCATION FOR IMPLICIT INTEGRATION. 

- LIMIT ON FREESTREAM EDDY VISCOSITY COMP. IN DFCFBL . 

- <UINF**2) 

- TIME (MINUTES) OF CPU USED. 

- (1.04) STEP SIZE MULTIPLIER. 

- STARTING TIME (DISTANCE). 

- (533.0) AIR REFERENCE TEMP. FOR COMPUTATIONS IN DIMEN 

- (533.0) REFERENCE TEMPERATURE. 

- (520.0) H2 REF. TEMPERATURE FOR COMPUTATIONS IN DIMEN 

- 1.0 + ( GAMMAF-1 . 0 ) * XMACHS**2 / 2.0 

- INITIAL LEVEL FOR SCALE FACTOR ' TSCALE ' 

- SCALE FACTOR FOR T2FIX * T3FIX AND T2PFIX IN PPRES. 

- STATIC TEMPERATURE COMPUTED IN CPINIT. 

- (12.0) LENGTH SCALE USED IN BRDSHW. 

- (2.0*PI> 

- ADD U2 CONVECTION TO RHS OF PP EQUATION. 

- ADD U2 * AND U3' TO RHS OF PP EQUATION. 

- ADD U3 CONVECTION TO RHS OF PP EQUATION. 

MASS WEIGHTED AVERAGE VELOCITY. 
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A 


RARRAY NAME DESCRIPTION 

ENTRY 

385 UCMULT - (1.0) CONVECTION TERM MULTIPLIER FOR U1 f TKEf EPS. 
203 UED - EDGE VELOCITY. 

63 UEDGE - UINF / UINFX USED IN BRDSHU. 

27 UINF - FREESTREAM VELOCITY. 

354 UINFX - (UINF) 

202 UWALL - VELOCITY JUST OF WALL. 

458 U1MIN - MINIMUM LEVEL OF U1 ALLOWED WHEN 'SETUP * IS CALLED. 

282 U2STRS - ADD REYNOLD STRESSES TO U2f U3 EQUATIONS. 

120 VARB - PACKED WORD OF VAR. BEING INTEGRATED. 

E.G. 12356.0 = U1 t U2f U3f Kf E 
348 VCMULT - CONVECTION TERM ADDED TO U2f U3 EQUATIONS. 

104 VELCST - UINF**2 / ( 2 . 0*G*AJ*CPA*T0A ) 

177 VLBTON - (1.488) VISCOSITY BRITISH TO MKS* 

178 VLBTOP - (14.88) VISCOSITY BRITISH TO CGS. 

346 VLDMLT - LAMINAR DIFFUSION ADDED TO U2f U3 EQUATIONS. 

102 VSTART - (101.0) PERCENT OF TD AT WHICH TO START 

U2f U3 COMPUTATIONS IN CONTES. 

383 WSMAX - WIDTH OF SOLUTION DOMAIN AT INITIAL STATION. 

273 XDELTA - (0.01) FOR EPSDIM .GT. EPSMAX * XDELTAf 

EPSDIM = EPSMAX * XDELTA. 

125 XLAM - (0.09) CONSTANT USED IN DFCFBL » 

109 XMA - (28.97) MOLECULAR WEIGHT OF AIR. 

61 XMACHO - MACH NUMBER. 

154 XMACHS - LOCAL MACH NUMBER. 

373 XMDOTC - AVERAGE MAS FLOW. 

66 XMF - (29.4) MOLECULAR WEIGHT OF FLUID. 

172 XMFACT - (UINF*SQRT ( XMA/ ( TOFINF*GAMMAF*G*RUNI V ) ) ) 

110 XMH - (2.016) MOLECULAR WEIGHT OF HYDROGEN. 

38 XMUINF - FREESTREAM VISCOSITY. 

98 XPRIME - NON-DIM. PRESUURE GRADIENT AT PRESENT STATION. 

52 XSCALE - (1.0) X1C0R SCALE FACTOR. 

11 XT - DIMENSIONAL STATION. 

201 XTC - PRESENT STATION FOR INTEGRAL PARAMETER PRINT. 

330 YMULT - (1.0) SCALE FACTOR FOR GRID MULTIPLIER. 

324 YNRMAD - (F10) DISPLACEMENT OF COORD. IN Y DIRECTION. 

198 YPLUS - Y+ VALUE AT WHICH TO SWITCH FROM MLT TO TKE ♦ 

53 YSCALE - (1.0) X2C0R SCALE FACTOR. 

382 ZMULT - (1.0) SCALE FACTOR FOR GRID MULTIPLIER. 

325 ZNRMAD - (G10) DISPLACEMENT OF COORD. IN Z DIRECTION. 
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